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Laser excitation of nanometer-sized atomic and molecular clusters offers various opportunities to 
explore and control ultrafast many-particle dynamics. Whereas weak laser fields allow the analysis 
of photoionization, excited-state relaxation, and structural modifications on these finite quantum 
systems, large-amplitude collective electron motion and Coulomb explosion can be induced with 
intense laser pulses. This review provides an overview of key phenomena arising from laser-cluster 
interactions with focus on nonlinear optical excitations and discusses the underlying processes 
according to the current understanding. A brief general survey covers basic cluster properties 
and excitation mechanisms relevant for laser-driven cluster dynamics. Then, after an excursion in 
theoretical and experimental methods, results for single- and multiphoton excitations are reviewed 
with emphasis on signatures from time- and angular resolved photoemission. A key issue of this 
review is the broad spectrum of phenomena arising from clusters exposed to strong fields, where 
the interaction with the laser pulse creates short-lived and dense nanoplasmas. The implications 
for technical developments include the controlled generation of ion, electron, and radiation pulses, 
as will be addressed along with corresponding examples. Finally, future prospects of laser-cluster 
research as well as experimental and theoretical challenges are discussed. 



PACS numbers: 36.40.-c,52.50.Jm 



Contents 



2. Above-threshold ionization and thermalization 



I. Introduction 

II. General survey of laser-cluster interactions 

A. Basic cluster properties and timescales 

B. Intense laser fields: key parameters 

C. Ionization and heating mechanisms in clusters 

D. Classification of coupling regimes 

III. Theoretical tools for cluster dynamics 

A. Approaches in general 

B. Effective microscopic theories 

1. The energy functional 

2. Time-dependent density-functional theory 

3. Semiclassical approaches 

4. Classical molecular dynamics 

C. Rate equations and the nanoplasma model 

IV. Experimental methods 

A. Generation of cluster beams 

B. Sources for intense radiation 

C. Particle detection techniques 

1. Determination of charge state distributions 

2. Acquisition of ion recoil energy spectra 

3. Energy and angular resolved electron detection 

V. Single- and multiphoton processes in clusters 

A. Single-photon electron emission 

1. Probing the density of states 

2. Angular distributions 

3. Time-resolved analysis 

B. Multiphoton signatures 

1. Competition of linear and nonlinear excitation 



VI. Cluster dynamics in strong fields 

A. Early surprises and basic trends 

1. Laser energy absorption 

2. Highly charged atomic ions 

3. Ion energy distributions 

4. Soft x-ray and EUV emission 

5. High harmonic generation 

B. Analyzing the microscopic cluster response 

1. The key role of collective excitations 

2. Difficulties of explaining high charge states 

3. Asymmetric ion and electron emission 

fpl VII. Perspectives of laser-cluster research 

A. Laser pulse shaping and control 

B. Towards VUV-, XUV, and soft x-ray pulses 

C. Clusters in an environment 

D. Relativistic particle acceleration with clusters 

E. Challenges for theory 



1£ 

m 
u 
u 

12 

y 

15 

m 

m 
m 

m 
m 

H 

19 

-1 I. INTRODUCTION 

m 
m 
m 
u 

1 

24 



Acknowledgments 



References 



Clusters of atoms and molecules frequently appear as 
a novel state of matter on the one-nanometer scale. For 
example, different types of bonding or various structural 
and chemical features can be realized within the same 
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material by just changing the particle size. The opportu- 
nity to vary, almost at will, the number of atoms in the 
clusters thus offers a unique avenue to explore the organi- 
zation and properties of matter from a fundamental point 



turbative nature. As a surprising featur e, clusters very 



efficiently absorb intense laser radiation (Ditmire et al. 



1997a), with an energy capture p er atom much highe r 



of view (Alonsc, 2006 


iHabcrland 




1994 




Martin 


7 


1996: 


Sugano and Koizumil 


1998). This also 


applies to opti- 



cal phenomena ar ising from small particles, e.g., due to 



surface plasmons ( Kreibig and VoUmcr 



cinated scientists since a long time (Mic, 1908; Raylcigh 



19951) which fas- 



1899). Today's lasers open an even more exciting per- 
spective of cluster science, i.e., the opportunity to steer 
and resolve ultrafast dynamics on the nanoscale. 
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Due t o the progress in laser technology (Keller 

2005), well-controlled short and intense laser 
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pulses can be routinely delivered these days. This opens 
the door to explore light-induced dynamical phenomena 
far beyond the mere analysis of ground state propcr- 



than for atoms or bulk material ( Batani et al. , 2001) 
Moreover, strong-field laser-clu ster interactions l ead t o 
the emission of fast ele ctrons (Bpringate et al. , 2003 ), 
multiply charged ions (pCoUcr et aI ^ 1999D , and high- 
energy photons ( McPhcrson et al. , 1994| ), documenting 
the excitation of core electrons. When compared to 
atoms, the appearance intensities for these products are 
strongly reduced with clusters. The discussion of the 
underlying dynamics and appropriate theoretical treat- 
ments is the central topic of this contribution. 

Different aspects of laser-excited clusters have previ- 
ously been reviewed, such as the electronic structure 



of simp l e me tal clusters (Brack, 1993; Ekardt 



dc Hcci , 1993), low- an d moderate-field dynamics (Rein- 



ties. For example, the real-time analysis of nuclear and hard and Suraud, 2003), ionizat ion mechanisms in stron g 



even electron motion becomes possible, as in the case of 
molecules or atoms (Corkuni and Krausz , 2007 ; Zewai^, 



1999 



199' ). When applied to clusters, short pulses controlled 
in amplitude and phase allow one to drive and resolve ion 
and electron dynamics on their natural time scales and 
under extreme conditions. For instance, electronic relax- 
ation processes or the time-evolution of collective modes 
can be studied with laser-excited clusters. As a more vi- 
olent scenario, strong-field exposure transforms clusters 
into well-isolated nanometer-sized plasmas, with inter- 
esting prospects for pulsed particle, radiation, or even 
neutron sources. With the advent of VUV free electron 



lasers (Feldhaus et al., 2005) coherent multiphoton inner- 
shell excitations are accessible with intense femtosecond 
pulses. Inspired by such opportunities, the subject of 
laser-cluster interactions has spawned sustained interdis- 
ciplinary activities and experienced enormous develop- 
ments over the last two decades. It definitely holds the 
promise to deliver unprecedented insights into the na- 
ture of light-matter interactions in complex systems and 
stimulated challenging efforts in experiment and theory. 

In this review we focus on the nonlinear response be- 
havior of clusters subject to laser fields, concentrating on 
the nonrelativistic intensity regime. It is our aim, in close 
connection between theory and experiment, to discuss 
signatures and mechanisms for multiphoton as well as for 
strong-field excitations. Nevertheless, even single-photon 
absorption can lead to complex dynamics, e.g., due to 
electron correlations, structural transitions, or compet- 
ing electronic decay channels. As a result, the response 
can clearly go beyond a simple and direct mapping of 
ground state properties. In any case, pronounced nonlin- 
earities emerge when multiphoton absorption is involved. 
As a typical example within the still photon-dominated 
regime, above-threshold ionization can be observed with 
clusters, showing additional finite-size and many-particle 
effects when compared to atomic systems. At higher in- 
tensities in the so-called field-dominated regime, the im- 
mediate excitation of several electrons and laser-driven 
collisions induce avalanche processes of highly nonper- 



optical and VUV laser fields ( Saalmann et al. , 2006 ) 
and exci t ations with ultraintense pulses (Krainov and 



Smirnov, 2002). The current report aims to deliver a 
present-day view on cluster dynamics in optical laser 
fields, with emphasis on the strong-field regime, and 
incorporates recent findings regarding angular resolved 
emission, electron acceleration, and processes behind 
very highly charged ions. Moreover, routes will be re- 
viewed to resolve the cluster response in time by varying 
the pulse duration or using dual-pulse excitations. Spe- 
cial features of this review are the extensive presentation 
of experimental and theoretical methods and the attempt 
of closely combining theory and experiment. 

The text is organized in six major parts. Section |l| of- 
fers a quick outlook of the topic and discusses basic phys- 
ical mechanisms. It thus provides some basic elementary 
stepping stones o n w hich to build an understanding of 
the topic. Section III is devoted to a brief survey of avail- 
able theoretical tools for describing cluster dynamics and 
tries to show how the various approaches may be linked 
together in terms of regimes for which they were primar- 
ily developed. Section IV focuses on experimental tech- 



niques, discussing cluster production and laser sources 
as a starter. Emphasis is essentially put on the detec- 
tion of emitted particles. In the ensuing presentation of 
selected results. Sec. ^ concentrates on the intermediate 
intensity domain in which photons still count. In this 
regime experiments have revealed detailed insight into 
the quantum nature of clusters and allow one to explore 
the emergence of nonlinear behaviors. Section VI finally 
comes to the main topic of the paper and describes highly 
nonlinear strong-field induced dynamics where quantum 
effects are partially wiped out. After a brief survey of 
initial/original results in the field, a detailed analysis of 
systematic trends and present day more elaborate ap- 
proaches are presented. This in particular concerns dif- 
ferential cross -sec tions and time-resolvcd analyses. Fi- 
nally, Section VII provides a brief outlook and proposes 
a few promising future directions of research in the field. 
We discuss in particular the prospects of laser develop- 
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mcnts, either in terms of pulse shaping of today's sources 
or by considering forthcoming devices hke the projected 
XFEL lasers. We also comment on embedded and de- 
posited clusters, high-energy particle acceleration with 
clusters, and point out some future challenges for theory. 



II. GENERAL SURVEY OF LASER-CLUSTER 
INTERACTIONS 

Laser irradiation of clusters allows the investigation 
of a broad spectrum of dynamical processes, ranging 
from single-photon driven ionization to the strong-field 
induced explosion of a nanometer-scaled plasma. Irre- 
spective of the regime under consideration, the absence 
of dissipation into substrate material offers a clean analy- 
sis of reaction products, i.e., electrons, ions, cluster frag- 
ments, as well as photons. Depending on the cluster ma- 
terial and the chosen laser intensity, quite different prop- 
erties and response mechanisms can be probed, as will 
be discussed throughout this review. Exemplarily, Fig. |l| 
illustrates a few response channels and properties that 
may be analyzed and can be viewed as a rough guideline. 

As an example for electron emission in the single- 
photon regime. Fig. [^a shows an ultraviolet photoelec- 
tron spectroscopy (UPS) result on Au^q obtained with 
low intensity laser excitation. The photoelectron energy 
spectrum images the electronic structure, i.e., binding 
energies and spectral occupation densities of single elec- 
tron states, and contains comprehensive information on 
the system. The large band gap in Fig. |l|a, for exam- 
ple, reflects the high stability of the tetrahedral Au2o (0 
et a I, |003 d). Besides structure analysis, photoelectron 
spectroscopy (PES) is a powerful tool for monitoring ex- 
cited states and reactions dynamics, see Sec. V.A . 

Laser- induced fragmentation may be analyzed, e.g., to 
determine optical properties. Fig. |l|b displays the opti- 
cal absorption cross-section of size-selected silver clus- 



kine tic ener gies by Coulomb explosion, see Sees. VLA.2 
and IVI.A.S . Examples for ion energy spectra from in- 
tense las er excitation of lea d clusters are displayed in 
Fig. |l|d ( [Peuber et ai , 2001 ) and document kinetic en- 



ters measured by photofragmentation ( Tiggesbaumker 
et c^, ^996, 1993). The spectra exhibit a pronou nced 
resonance, i.e., the Mie surface plasmon, see Sees. ILA 



and |LC. Collective excitations, as prime examples for 
multieleetron effects, are not only relevant in the single- 
photon limit, but are important for the cluster response 
in th e mult ipho ton a nd strong-field regime as well, see 
Sees. IV.B.ll and IVLB. 



With increasing laser intensity, nonlinear and feedback 
effects begin to severely influence the cluster response, 
such as the electron emission. Fig. |l|c shows an example 
for larger silver clusters, where the measured total elec- 
tron yield, i.e., the average cluster ionization, is plotted 
as a fu nction of the te mporal width of the exciting laser 
pulse ( Radcliffc , 2004). The strong variation with pulse 
duration reveals a pronounced ionization dynamics that 
can be related to the interplay of collective plasma heat- 
ing and ult rafast relaxation of the ionic structure, see 
Sec. VLB.l , In addition, as a result of high charging of 
cluster constituents, atomic ions are accelerated to high 



ergies of up to hundreds of keV as well as a clear clus- 
ter size effect in the recoil energy. Within the strong- 
field induced excitation process a hot and highly ionized 
nanoplasma is formed. Clear evidence for the presence 
of energetic electrons is given by the creation of inner- 
shell atomic vacancies in the cluster constituents, the 
recombination of which can be monitored by analyzing 
the extreme ultraviolet (EUV) and x-ray emission, see 
Sees. VI. A. 4 and VI.B.2| . The example in Fig. |l|e shows 
energetic Xc(-radiation at 12.6 keV resulting from irradi- 
ation of krypton clusters ( Issac et ai , 2004). A detailed 
analysis of the EUV and x-ray emission can be used for 
monitoring ion charge state distributions. 

The examples highlighted in Fig. |^ illustrate the wide 
spectrum of phenomena resulting from laser irradiation 
of clusters. Before analyzing particular response effects 
in more detail, a few basic facts about "protagonists" of 
such processes, i.e., clusters and lasers, will be recalled. 
In the following we furthermore remind basic mechanisms 
of energy absorption and ionization relying on both in- 
dividual atomic and cooperative processes and provide a 
rough classification of different coupling regimes. 



A. Basic cluster properties and timescales 

Cluster properties are strongly dependent on the type 
of their constituents. We consider four typical cluster 
materials: Na as a simple metal, Ag as a noble metal, 
C as a eovalent material, and Ar as a rare-gas systems. 
Table | recalls a few basic facts of these elements, e.g., 
the electronic core and valence levels and corresponding 
energy gaps. Since cluster properties are by nature also 
size-dependent (number of constituents between a few 
and several thousand atoms), atomic, dimer, and bulk 
values are stated, which fixes typical orders of magnitude. 

For a given element, the atomic ionization potential 
(IP) and the bulk work function (WF) indicate the elec- 
tronic stability of a corresponding atomic cluster with 
respect to optical excitation. Both IP and WF follow a 
similar trend over the given materials, i.e., increase from 
Na to Ar. Typically, metal clusters can be ionized or ex- 
cited much easier, i.e., with lower photon energies or less 
intense radiation, than eovalent or rare gas systems. This 
trend is also reflected in the first atomic dipole transition 
(lowest dipole excitation). The IP further indicates the 
ionization behavior in strong fields as it determines the 
critical laser inte nsity required for atomic barrier sup- 
pression, see Sec. [I.C for details. 

Structural stability is not necessarily linked to that of 
the electronic system. This becomes evident after com- 
paring dimer dissociation energies or bulk cohesive en- 
ergies with the IP's, e.g., for C with Ar. Note that the 
bulk cohesive energies roughly reflect the binding energy 
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FIG. 1 Five decay channels of laser-excited clusters aside with properties/processes that may be resolved from their analysis 
(see text), (a) electr onic structure o f negatively charged gold clusters with 20 atoms (Aujq) extracted from the photoelec- 
tron spectrum, from ( Li et al , 20Q3a|), with permission i: r om A A AS; (b) optical absorption of Ag^ and Ag^ as determined by 



photofragmentation, adapted from ( Piggcsbaumker et al., 1996, 1993), with permission from Elsevier; (c) ionization dynam ics of 



Agjv i n int ense laser pulses resolved by measuring the total electron yield as function of pulse width at fixed pulse energy (iRad 



cliff c| 2004| ); (d) Coulomb explosion of Pb]v analyzed by recoil energy spectroscopy of emitted atomic ions, from (Teuber et al. 



200 1|), with kind permission of The European P hysical Journal (E PJ); (e) inner-shell recombination in strongly excited krypt' 



clusters measured by x-ray spectroscopy, from (Issac et al, 2004), with permission from American Institute of Physics. 





Na 


Ag 


C 


Ar 


atom 




ionization potential'^' 


[eV] 


5.14 


7.58 


11.26 


15.8 


e , - £ 111 


[eV] 


26.0 


53.9 


8.21 




valence level 




3s 


5s 


2p 




core level 




2p 


4p 


2s 


3p 


lowest dipole exc'^' 


[eV] 


2.1 


3.66 


7.48 


11.62 


critical laser intensity 


L cm- J 


3x10^^ 


1x10" 


6x10" 


2x10" 


dimer 




bond length'2'3'*'^' 


[A] 


3.08 


2.53 


1.20 


3.83 


dissoc. energy'^'^' 


[eV] 


0.76 


1.69 


6.3 


0.012 


bulk 




work function'^' 


[eV] 


2.75 


4.26 


4.8 


15.8 


cohesive energy'^' 


[eV] 


1.12 


2.95 


7.8 


0.08 


Wigner-Seitz radius'^' 


[A] 


2.10 


1.59 


1.21 


2.21 



TABLE I Basic atom, dimer, and bulk properties for four 
typical cluster materials. Bulk properties for C correspond 
to graphite which is close to the Ceo cluster and carbon nan- 
otubes. The c ritical laser intensity is estimated with Eq. 
see Sec. [I.C. The Wigner-Se itz rad i us r, ch aracterizes the 
atom ic density. '^N I ST; P '( |Weast|, [198^);'^' (|Verma et al 
198^);'*i(|Be"utel et al\, |l993|);'^' ([Hirsclifelder et o/.|, |l954[). 



per atom of the cluster, while the atomic Wigner-Seitz 
radius Tg may be used to approximate the cluster radius 
{Rci ~ TgiV^/'^). The values for the dimer bond length 
indicate typical interatomic distances. 

In the visible and ultraviolet spectral range the optical 
response is mainly determined by valence electrons. In 
metal clusters, electron delocalization leads to a strong 
resonance, the Mie surface-plasmon, as a unique feature 
of finite objects with sub-wavelength dimension. It cor- 
responds to a collective oscillation of the whole valence 
electron cloud against the ionic background. When con- 



sidering schematically a cluster as a metallic drop ( Mie , 
1908| ), the Mie surface plas mon frequency of a neutra l 



system can be estimated as ( Brack) , [1993 ; de Heei, 1992) 



e (47reomer^ 



-1/2 



(1) 



with rs the effective Wigner-Seitz radius of conduc- 
tion electrons, e the elementary charge, eo the per- 
mittivity of vacuum, and mg the electron mass. For 
small NaAT, for example, the plasmon en er gy is around 
fiujMic ~ 2.8eV ( ^chmidt and Haberland , 1999), while 
Eq. (0) predicts a value of 3.4eV. This indicates that 
the actual Mie response depends on further details (finite 
size effects, geometrical structure, excitation, net charge, 
etc.), but Eq. (|l|) already provides a reasonable order of 
magnitude sufficient for many forthcoming discussions. 
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electron-electron / electron-ion collisions 
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electron-electron collision time: Tcc ~2-3fseV l{kBT) 
electron evaporation time: ^cvap~ {^bTY 

FIG. 2 Typical time scales for the dynamics, taking sodium 
clusters as a prototype. On the top the ranges associated to 
fs lasers are depicted . Further, processes related to motion 
(cycle times) and lifetimes due to relaxation (decay times) are 
indicated. Approximate expressions for electron-electron col- 
lisions (tcc) and electron evaporation (rcvap) are given below. 



For considering reaction pathways and energy dissipa- 
tion it is useful to compare relevant time scales. To that 
end we consider Na as a typical example for a metal clus- 
ter. Fig. ^ provides a schematic overview over times re- 
lated to laser characteristics, electronic and ionic motion, 
and lifetimes for relaxation processes. For the moment 
we ignore the extremely short times associated with core 
electrons. They certainly play an important role in in- 
tense laser fields, but are usually dealt with in terms of 
simplified rate equations, see, e.g., Sec. [II. C . The pulse 
duration of optical lasers may be varied over a wide range 
extending from fs to ps or even ns. We focus here on pulse 
widths of the order a few tens to a few hundred fs. 

The shortest time scales in Fig. || are related to the 
electronic motion. The Mie piasmon period as the most 
basic one is of the order of fs, cf. Eq. ([^). In the 
same range, but with a wider span from sub-fs to sev- 
eral fs, are cycle times for other single-particle excita- 
tions and direct electron escape, i.e., single-particle ex- 
citation into the continuum. Somewhat slower is the 
piasmon decay due to Landau fragmentation, in analogy 



to Landau damping k nown from plasma physics (Lifs- 
chit:; and Pitajewski, 1988). In clusters. Landau frag- 



mentation results from the coupling of plasmons with 
energetically close single-particle excitations. Viewed 
in coordinate space, it corresponds to collisions of elec- 
trons with the anharmonic potential at the cluster sur- 
face. The Landau relaxation time depends on clus- 
ter size and has, e.g., for Najy, its lowest values for 



N « 1000 dBabst and Reinhar4 |1997| ). For N > 1000 it 
can be estimated from the time between collisions of an 
electron with the cluster boundary ("wall friction") as 
« {rsN^^^)/ v^ where v„ = (n / m)(97 r/4)^/^/rs is the 
Fermi velocity ( [Yannoulcas et al , [1990D . For N < 1000, 
however, increases for smaller N due to the reduced 
level density. The relaxation time describes damping 
due to electron- ion collisions. It is strongly temperature 
dependent (~ 30 fs for Na at 273K) and scales as Tei oc 



at low temperature due to e lcctron-phonon scatter- 
ing ( Ashcroft and Mermin , 197 1 ) and f ollow s oc T^/^ 
in a high-temperature plasma (Spitzer, 1956). 

The most widely varying times are related to the colli- 
sional damping from electron-electron collisions and ther- 
mal electron evaporation. Both strongly depend on the 
internal excitation of the cluster, which may be charac- 
terized by an electronic temperature T. A simple con- 
nection between internal excitation energy per electron 
e* = E* /N and temperature can be estabhshed by the 
Fermi gas model. For k^T <^ e^, T can be estimated as 
kBT = 2{e^e*y/^/n, where = {971^4)^/'^ /{2myj 
is the Fermi energy. For the particular case of sodium at 
bulk density, we have ksT ^ (1.28 eV e*)^/^ Electron- 
electron collisions are the key mechanism for electronic 
thcrmalization. The law for the corresponding col- 
lision tnn£_Jii_^j£^_^|_2S_faiowi^ 

ory ( Kadanoff and Baym, 1962; Pines and Nozieres, 
1966| ). For low T, collisions are strongly suppressed 



due to Pauli blocking of energetically available electronic 
states. At high T, electron collisions become competitive 
with Landau damping and sometimes even the dominat- 
ing damping mechanism. Electron-electron collisions can 



be described by semiclassical models, see Sec. III.B.3 



An even more dramatic temperature (or excitation 
energy) dependence appears for the electron evapo- 
ration time, whose trend is dominated by the expo- 
nential factor exp {Eip /(kBT)), where Eip denotes the 
value of the ionization potential. The more detailed 
expression for the evaporation time given in Fig. |2| 
is based on the Weisskopf formulae ( Weisskop^ 1937| ) 
Tevap « ^?iV(8"ier2iV2/3)(fcBT)-2 exp {Eip/{kBT)) and 
a cluster size of A^=100. For this size the crossing 
point Tcou ~ Tevap occurs at a temperature of about 
ksT = 0.8 eV. This corresponds to a hot (nano-)plasma 
where finite electron clouds are practically an unstable 
evaporative ensemble. In general, electron evaporation 
represents a very (sometimes even the most) efficient 
cooling mechanism for highly excited clusters. 

Ionic motion spans a wide range of long time scales. Vi- 
brations, which may be measured by Raman scattering, 
see, e.g., (Portales et al., 2001), are typically in the meV 
regime, i.e., have cycle times of 100 fs to Ips. In small 
clusters, ionic v ibrations can induce satellites in t he op - 
tical spectrum (EUcrt et al, ^002 ; Fehrer et al, 2006). 
Strong laser irradiation usually leads to large ampli- 
tude ionic motion and cluster explosion due to Coulomb 
pressure generated by ionization and thermal excitation. 
Electron-ion coupling due to Coulomb pressure proceeds 
at the electronic time scale, i.e., within a few fs. The 
effect on the ions, however, develops at slower scale, typ- 
ically beyond 100 fs, due to the large ionic mass. The 
time scale of Coulomb explosion can be estimated by 
considering sudden ionization of cluster constituents to 
an average atomic charge state (g) . In this case the clus- 
ter expands homogenously and doubles its radius after 
Tdoub ~ 2.3(V27reo/e) m\ilrV^ / (q), where mion is the ion 
mass and is the initial atomic Wigner-Seitz radius. For 
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NaAT this yields Tdoub ~ 63fs/{q). In consequence, strong 
ionization drives clusters apart quite rapidly, accompa- 
nied with strong changes in the optical properties. Cor- 
responding signature s can be analyzed with pump-probe 
techniques, see Sec 



Field strength [V/A] 

10 100 1000 



VLB 



For excitations that do not 
induce explosion, the time scale of electron-ion thcrmal - 
ization reaches up to the ns range ( iFehrcr et al\ , ^006\ ). 
Ionic relaxation is even slower, e.g., thermal emission of 
a monomer can easily last /iS. 

As shown above, cluster dynamics comprises a large 
span of time scales, making their theoretical description 
to a great challenge. Ionic motion may require a simu- 
lation time up to several ps while electronic times scales 
down to a small fraction of a fs have to be resolved. The- 
oretical approa che s for a corresponding description are 
subject of Sec. III. Relaxation processes at the ns scale. 



however, require more phcnomcnological approaches. 



B. Intense laser fields: key parameters 

We proceed with a brief summary of basic facts and 
key parameters of intense laser fields. In the nonrela- 
tivistic regime, laser pulses acting on atoms, molecules, 
or clusters can usually be described as a homogenous 
time-dependent electric field of the form 



\£o\f{t) cos(a;iast + 



(2) 



where denotes linear polarization in z-direction, £q 
is the peak field strength, f{t) is the normalized tem- 
poral field envelope of the pulse, fkJus is the photon 
energy of the carrier, and (p{t) is an additional tempo- 
ral phase. Any other polarization (linear or circular) 
can be described by superposition. The phase can be 
written as (p{t) = ipce + f + 3*'^ + 0{t'^) where (pcc 
is the carrier-envelope phase, P and 7 denote linear and 
quadratic chirp, and the last term indicates higher or- 
der chirp contributions. Furthermore, the instantaneous 
frequency reads WinstlO ~ ^las + ^(t) and the instanta- 
neous pulse intensity is given by I{t) = /o|/(f)p, where 
Iq = ceo|5oP/2 is the peak intensity and c is the vac- 
uum speed of light. Typically, the pulse duration r is 
given as the full width at half maximum (FWHM) of the 
temporal intensity profile. A common temporal pulse 
profile is a Gaussian field envelope, which then reads 
f{t) = exp (— 21n2t^/T^). In absence of chirp, the band- 
width Au! (FWHM) of the corresponding spectral inten- 
sity profile is related to the temporal pulse width via the 
time-bandwidth product toAuj/2t: = 0.441. Increasing 
the pulse duration by dispersive pulse stretching to r in- 
duces a linear chirp of /? = ±41n2\/s^ — 1/(s^Tq), where 
s ~ t/tq > 1 is the stretching factor with respect to 
the bandwidth-limited pulse. The chirp direction (up or 
down) depends on the sign of the group velocity disper- 
sion of the optical element. However, it should be noted 
that the exact forms of f{t) and (p{t) are not always easy 
to ascertain experimentally. Nonetheless, the pulse du- 
ration can nowadays be varied very flexibly over a wide 
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FIG. 3 Intensity-frequency regimes attainable with different 
liigh intensity laser systems (shaded blocks). Corresponding 
wavelengths and electric field strengths are displayed on the 
additional scales. Lines indicate regions of constant pondero- 
motive potential Up. The transition from photon- to field- 
dominated coupling is roughly given by Up = Eip , as schemat- 
ically depicted for an IP of a few eV. VUV-FEL: vacuum ul- 
traviolet free electron laser; X-FEL: x-ray free electron laser. 



range, e.g., between a few fs up to ns for optical lasers. 

In the dipole approximation and using the length 
gauge, the coupling of the pulse to an electron at po- 
sition r can be described by an external potential 



(3) 



Therefore the system size has to be well below the wave- 
length A ~ 27rc/a;ias, which is well justified for nm clus- 
ters and excitation in the optical domain (A ~ /xm). The 
dipole approximation becomes questionable for UV pho- 
tons and very large clusters, but will be valid in most 
cases considered below. 

To classify coupling regimes it is useful to consider a 
freely oscillating electron (pure quiver motion, no drift 
velocity) in the laser field. The cycle averaged kinetic 
energy defines the ponderomotive potential, which reads 



(4) 



at the pulse peak. It can be expressed more conve- 
niently by Up = 9.33 X 10-"eV x /o[W/cm2] (A[^m])2. 
Fig. 1^ displays the dependence of Up in the frequency- 
intensity plane aside with the characteristic parameter 
regions which can be realized with high intensity laser 
sources. As a rule of thumb, regimes of photon- and field- 
dominated coupling are separated by a Up that equals the 
typical electron binding energy in the considered system, 
as schematically displayed in Fig. ||. This condition is re- 
lated to the Kcldysh parameter, as discussed in more de- 
tail in Sec. [I.C. Figure ^ further illustrates the enormous 
flexibility of optical lasers to produce high intensities up 
to the relativistic limit where Up becomes nonnegligible 
compared to the electron rest energy. In this review, 
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however, we focus on intensities for which relativistic ef- 
fects and the magnetic field of the pulses may be ne- 
glected. Compared to optical lasers, vacuum ultraviolet 
and x-ray free electron lasers (VUV-FEL/X-FEL) cover 
a fundamentally different regime, i.e., photon-driven dy- 
namics at high intensities due t o the low ponderomotive 
potential, sec ( ^aalmann et al , 2006| ) and Sec. VII. B, 



C. Ionization and heating mechanisms in clusters 

Several basic ionization and energy absorption mecha- 
nisms arc of relevance for describing laser irradiated par- 
ticles and will be briefly introduced below. Departing 
from concepts for atomic and molecular systems wc move 
on to cooperative and collective effects which stem from 
the many-particle nature of clusters. 

On the atomic level, two fundamentally different pho- 
toionization processes may be considered. The first is 
vertical excitation of a bound electron by single- or mul- 
tiphoton absorption in a rapidly oscillating laser field, 
see multiphoton ionization (MPI) in Fig. This mech- 
anism proceeds over many laser cycles and prevails for 
weak and moderate fields in the so-called perturbativc 
domain. A MPI process of order i/ is characterized by 
the reaction rate T^, = <t,jI'^ , where Ui, is the correspond- 
ing cross-section. MPI, which may be enhanced when 
intermediate resonant states are available, can promote 
electrons far beyond the continuum threshold, leading 
to characteristic peaks separated by units of the pho- 
ton energy in the electron energy spectrum. This ef- 
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FIG. 4 Schematic view of ionization mechanisms in atoms 
molecules and clusters. Panels (a) and (b) display potentials 
of the unperturbed ions Vlon, the laser Vias, and their effective 
sum. In panel (a) the pathways for MPI and OFI of a bound 
electron are indicated, while panel (b) depicts the charge- 
resonance-enhanced ionization (CREI). The vertical arrows 
in (b) indicate the Stark shift. Panel (c) illustrates inner and 
oute r ionization of a cluster based on an effective potential , — 



feet, termed as above-threshold ionization (ATI), is well- 
kno wn from atoms and also appears in clusters, see 
Sec. V.B.2| . The second mechanism is optical field ion- 
ization (OFI). Here the laser acts as a quasistationary 
electric field. For sufRciently strong fields, bound elec- 
trons tunnel through the barrier emerging from the com- 
bined potential of the residual g-charged ion and the laser 
field, i.e., V{x) oc — a/jzl — z, with a = (ge^)/(47reo|^-o|)- 
This is schematically depicted in Fig. ^ (dashed curve). 
The probability for atomic tunneling ionization can be 
described by the well-known A DK rates found by Am - 
mosov and Delone and Krainov (Ammosov et al, 1986). 

A useful measure for the significance of MPI over O FI 
is the Keldysh adiabaticity parameter ( Keldysh , 1965| ) 
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/ Eip 



(5) 



which compares the IP with the peak kinetic energy of 
a freely quivering electron (2C/p). Single- or multiphoton 
ionization dominates for 7 3> 1, where the quiver energy 
is small compared to the IP. For 7^1, the binding en- 
ergy can be overcome within a single laser cycle and OFI 
is promoted. An equivalent expression for the Keldysh 
parameter is 7 = wiasTtunnei, which gives a ratio of the 
tunneling time Ttunnci = ■\/2i?ip nie/ {c'^Eq) and the opti- 
cal period. Optical field ionization dominates if the tun- 
neling time is comparable to or smaller than the optical 
period; MPI is the leading process otherwise. 

Within the tunneling regime (7 ^ 1), the ioniza- 
tion probability in one optical cycle approaches unity if 
the potential barrier can be fully suppressed. For an 
atomic system, this so-called barrier suppression ioniza- 
tion (BSI) roughly sets in at the threshold intensity 



2e6 



o%^4xlO^ (^^P[f]) W/cm^ (6) 



which reason ably predicts ion a ppearance intensi_ties in 
atomic gases (Angst et al., 1989). Note that Eq. 



was 



used to determine the critical intensities in Tab 

The above considerations apply to isolated atoms 
where the laser parameters govern the dynamics. For 
extended systems, i.e., from the molecular level on, struc- 
tural details become increasingly important. Ionization 
barriers are influenced by the fields from neighboring 
ions, which, for example, gives rise to charge-resonance- 
enhanced ionization (CREI) well- known from strong-field 
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ionization of diatomic mo lecules (Seidcman et al. 
Zuo and Bandrauk , 199E ). Within this process, an ap- 
propriate internuclear separation results in a simultane- 
ous lowering or suppression of inner and outer potential 
barriers with respect to the Stark-shifted electronic states 
(see Fig. ^), giving rise to an enhanced ionization rate. 
For larger or smaller separations either the inner or the 
outer barriers increase and the ionization probability is 
reduced. As a truly cooperative effect, CREI has been 



considered also for very small clusters (Sicdschlag and 



Ros^ , |200^ ; [Veniard et al] , |200T]), cf. Scc. |VI.B.l 
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Very convenient for describing charging dynamics in 
larger systems is the concept of inner and outer ioniza- 
tion (Last and Jortncr, 1999). As indicated in Fig. 
electrons in the cluster may be classified into tightly 
bound, quasifree, and continuum electrons. Within 
this picture, inner ionization describes the excitation of 
tightly bound electrons to the conduction band, i.e., elec- 
trons are removed from their host ion but reside within 
the cluster. Correspondingly, the final excitation into the 
continuum is termed outer ionization, which contributes 
to the net ionization of the system. At moderate laser in- 
tensities, systems with initially delocalized electrons, like 
metallic particles, may undergo outer ionization only. In 
any case, however, the energy span between the thresh- 
olds for inner and outer ionization grows with cluster 
charge, cf. Fig. underlining the growing importance 
of quasifree electrons for the interaction dynamics. Be- 
sides purely laser-induced MPI and OFI, ionization can 
be driven by cluster polarization (field amplification) or 
cluster space-charge fields, e.g., subsequent to strong ion- 
ization. In addition, quasifree electrons can drive electron 
impact ionization (EII ) , as may be described by semiem- 
pirical cross-sections (Lotz, 1967). The onset and self- 
amplification of such addit ional processes is f r equen tly 
termed ionization ignition ( Rose-Petruck et ai, 1997). 

The presence of a nanoplasma, i.e., of quasifree elec- 
trons and (multi-)chargcd atomic ions in the cluster, 
has substantial impact on the energy capture from a 
laser pulse. If collective effects are negligible, elec- 
trons can acquire energy from the laser field via Inverse 
Bremsstrahlung (IBS), i.e., by absorbing radiation en- 
ergy during scattering in the Coulomb field of the ions. 
IBS relics on the conversion of laser-driven electron mo- 
tion into thermal energy because of directional momen- 
tum redistribution within elastic collisions and is a basic 



volum e-heating effect in underdense plasmas (Krainov, 
200C ) . Considering a fixed coUisional dcphasing time Tcoh 
(inverse collision frequency), the IBS heating rate per 
electron in terms of the ponderomotive potential reads 



dE\ ^2j,_W^ 



coU^las + 1 



(7) 



Whereas the heating rate becomes independent of wias in 
the low-frequency case (dc-limit), a C/p/rcoii-dependence 



is found for tr 



^1. It should be noted that the 



coUisional relaxation time, which is a function of electron 



temperature (cf. Sec. [LA) and becomes frequency de- 
pendent (tcou oc uj^l^) for short- wavelength laser excita- 
tion, is in general difficult to obtain. For laser-irradiated 
clusters, pure IBS heating dominates the energy capture 
of quasifree electrons only at laser frequency far above the 
Mie plasmon frequency. If the laser frequency becomes 
comparable to or smaller than WMic, the collective re- 
sponse of quasifree electrons in the cluster has to be taken 
into account. Surface charges from the laser-driven col- 
lective electron displacement induce polarization fields, 
that strongly modify the effective field in the cluster in 
amplitude and phase. For a spherical plasma and suffi- 



ciently small displacements the corresponding restoring 
force is linear, i.e., the absorption rate per electron for 
collective IBS heating is described by a Lorentz profile 



dE\ 

It /„ ~ 

I Res 



Tcollt^L 



'coll ' 



/ ? 9 \ ^ , 9 



-Tas 



(8) 



This expression is equivalent to the heati ng rat e assumed 
in Ditmire's nanoplasma model, cf. Sec. JIII.C . Whereas 
the absorption rates in Eqs. (0) and (g) meet in the 
high-frequency limit, IBS heating is strongly suppressed 
for wias <C wmio due to efhcient screening of the external 
field by the collective electron displacement. Most im- 
portantly, excitation with lo\^ k, wmic leads to plasmon- 
enhanced energy absorption in Eq. (^), cf. the cross- 
sections in Fig. |l|b. Resonant collective driving of clus- 
ter electrons can produce strong field amplification that 
supp orts cluster ionization and direct acceleration ofclcc- 
tron ( Fennel et at , 2007a ; Rcinhard and Suraud , 199^ ). 

In the above discussion the absorption rates have been 
assumed to scale linearly with intensity (cx Up), cf. 
Eqs. (|^) and (^). This requires that the dephasing time 
and the plasmon frequency are constants. In strong 
fields, however, the large quiver amplitudes actively mod- 
ify the nanoplasma properties. Hence, both the dephas- 
ing time and the plasmon frequency become functions of 
intensity which introduces additional nonlinear terms. 

Another very important aspect for the cluster response 
to strong optical laser fields is the time dependence of 
the plasmon energy. It scales as wmic oc ^/phg{t), where 
Pbg is the ion-background charge density. In early stages 
of the interaction pbg is usually too high for being in 
resonance with the driving IR-field, i.e., the system is 
overcritical. This is the case in metal- and, already after 
moderate inner ionization, in rare-gas clusters and leads 
to strongly suppressed IBS heating as explained above. 
Less efficient surfa ce heating effects like vacuum hea ting 
or Brunei-heating (Brunei, 1987; Taguchi et al, 2004) re- 
main active in this overcritical state. Therefore electrons 
that are pulled away from the surface by the laser field are 
accelerated outside and contribute their acquired energy 
upon recoUision with the cluster. In any case, as a result 
of moderate charging and heating. Coulomb forces and 
thermal electron pressure eventually induce an expansion 
of the cluster. Corresponding time scales a re ty pically 
between a few tens of fs to some ps, see Sec. 11. A for an 



estimate of the radius doubling time for pure Coulomb 
explosion. With cluster expansion the frequency of the 
collective mode decreases and transiently matches the 
laser frequency at a certain time, producing a short-lived 
but strong absorption enhancement, cf. Eq. (j|). This 
idea is a central eleme nt of the hydrodynamic approach 
from Ditmire, see Sec. III.C, however, characterizing the 



resonance condition in terms of a critical electron den- 
sity. The latter is justified only for nearly charge neutral 
systems, such as very large clusters. Since, according to 
the harmonic potential theorem ( pobson , 1994), the ionic 
background creates the restoring force for quasifree elec- 
trons, the background charge density is the more general 
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parameter applicable also to charged systems. Nonethe- 
less, for sufficiently long pulses the transient resonance 
induces efficient heating of quasifree electrons and, as a 
consequence, strongly supports outer ionization and clus- 
ter Coulomb explosion. At high laser intensity, this de- 
layed resonant coupling is important irrespective of the 
cluster material and leaves clear signatures i n the ab sorp- 
tion as well as in emission spectra, see Sec. VI.B.l . 
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D. Classification of coupling regimes 

While the relative importance of the above mechanisms 
depends on the specific scenario, regimes can be identified 
where particular processes prevail. However, such classi- 
fication cannot be achieved based on a single parameter 
like laser intensity. While very low intensities lead to lin- 
ear and very high ones to nonlinear behavior, other laser 
characteristics or cluster properties determine the nature 
of the response for intermediate cases. We briefly discuss 
a rough sorting of regimes used throughout this article. 

The linear regime is the domain of weak laser fields 
associated with single-photon processes and large values 
of 7 (cf . Eq. . Mechanism arc sensitive predominantly 
to the laser frequency. The prevailing examples are op- 
tical response spectra. As this is a key to ol, the r e is a 
huge body of rev i ews and b ooks s ee, e.g., ( Brack 1993[ 
Haberlaiiq , f994; ie Heer, f993; Kreibig and VoUmeij , 
1995| ). Early cluster experiments often used ns pulses 



for studies on structure or low -energy dynamics (Habcr- 
lan4 |1994| ; [Naher et al\ , |1997| ). Another typical process 



is single-photon ionization which can be analy zed b y pho- 
toelectron spectroscopy, see Fig. 



V.A 



and Sec. 

The muUiphoton regime is associated with moderate 
laser intensities where multiphoton processes begin to 
show up (/ ^ 10^ — 10^"^ W/cm^ depending on material 
and frequency). Each laser parameter, i.e., frequency, 
field strength, and pulse profile, becomes equally im- 
portan t. Typical e xamp l es are second harmonic g cner- 



ation ( iGotz et all |1995| ; [Klein- Wiele et alj |1999D and 



multiphoton ionization. Of particular interest are cases 
where a multiple of the photon energy can excite an inter- 
mediate state of the system. Then, besides direct MPI, a 
sequential ionization fro m the (long -living) intermediate 
state becomes possible (Pohl et at, 200f ). Another ex- 



ample is above-threshold ionization. Processes emerging 
in the multiphoton regime are subject of Sec. V.B. 

At sufficiently high intensity the laser irradia- 
tion produces large ionization and strong heating 
(/ ^ 10^^ — 10^^ W/cm^). The excitation of many elec- 
trons and strong feedback effects on the response indi- 
cate the so-called strong-field domain where the dynam- 
ics cannot be treated perturbatively. Typically, the ex- 
citation leads to cluster Coulomb explosion, accompa- 
nied by emission of energetic particles, i.e., electrons and 
ions, as well as photons. The emitted ions usually carry 
higher charges than in the case of irradiation of single 
atoms which underlines the impact of cooperative pro- 
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FIG. 5 Ionization of Nag" as function of laser intensity for 
excitation by 70 fs cos^-shaped laser pulses for two frequencies 
(as indicated). The ionization potential is 7.2 eV. Three pho- 
tons of Tiojias = 2.7 eV are required to lift an electron into the 
continuum (multiphoton ionization) while one photon suffices 
for ?!.Wias = 10 eV (linear behavior). At high intensity both 
cases become nonperturbative, indicating strong-field condi- 
tions. Note, that Tilji^b = 2.7 eV is close to the Mie plasmon 
of Na,g , which leads to the early onset of the strong-field re- 
sponse in this case. Calculations are done in TDLDA. 



cesses. Moreover, the reactions proceed somehow similar 
for very different cluster materials (from metals to rare 
gases) since electrons from atomic shells are activated and 
the transient nanoplasma determines the dynamics. Such 



highly nonlinear processes are in the focus of Sec. VI 



A possible marker for the actual regime is the total 
ionization yield as function of laser intensity. Lowest 
order perturbation theory predicts that the yield scales 
with cx where v is the number of photons required to 
overcome the ionization potential. Figure ^ gives an ex- 
ample for Nat}' excited with 70 fs laser pulses and shows 
the intensity-dependent electron yield for two different 
laser frequencies. The slope at low intensities agrees 
nicely with the I'^ law, yielding v ~ 3 (multiphoton) 
for the lower- and v — I (single-photon) for the higher 
frequency. However, the curves turn over at higher inten- 
sities where sorting in orders of photon becomes obsolete 
(breakdown of perturbation theory) . One approaches the 
"strong-field domain". Note, that the two laser frequen- 
cies perform in a very different way. With Ticjias = 10 eV 
excitation, the yield follows the linear behavior and be- 
comes nonperturbative at rather large intensities. With 
the lower frequency the ionization is a three-photon pro- 
cess and the transition to the nonlinear regime evolves 
at a much lower intensity. Two effects contribute in the 
latter case: the near-resonance exc itation of the Mie plas- 
mon (Reinhard and Suraud, 199S) and the stronger im- 
pact of optical field effects at lower Keldysh parameters. 
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TABLE II Hierarchy of approaches for the description of electrons and ions in a cluster. Acronyms are defined in the text. The 
range of applications is listed in the column regime where structure is abbreviated as S, excitation spectra (optical response) as 
E, and dynamics as D. The label D* indicates the capability to describe electron emission and E* stands for excitation energy. 



III. THEORETICAL TOOLS FOR CLUSTER DYNAMICS 
A. Approaches in general 

Clusters are complex systems and their theoretical de- 
scription requires approximations to the full quantum- 
mechanical many-body problem - the more so for truly 
dynamical situations. As approximations are always a 
compromise between feasibility and demands, there ex- 
ists a rich spectrum of methods. Table ^ tries to provide 
a brief overview of commonly applied methods - in the 
upper part for electrons and in the lower part for the 
ions. Keywords, numbers, and citations are guidelines 
and by no means exhaustive. They should be under- 
stood as examples and estimates of orders of magnitude. 
For ab-initio methods some entries for typical sizes and 
excitation energies E* are left open as they have, in prin- 
ciple, a huge range of validity, but are, in practice, very 



limi ;cd by quickly growing numerical expense. We add a 
few remarks while going through the table. 



T [ie class of ab initio theories covers a huge range of 
treatments depending on the size of the underlying ba- 
sis space, in particular for the configuration interaction 
(CI) and the multiconfiguration time-dependent Hartree- 
Fock (MC-TDHF) approach. The most general methods, 
i.e., exact time-dependent Schrodingcr equation (TDSE) 
and Quantum Monte-Carlo (QMC), are still restricted to 
very few electrons and presently not applicable to clus- 
ters. The vast majority of theoretical investigations of 
cluster dynamics with quantum aspects relies on density- 
functional theory (DFT) based methods, with quantum 



mechanical (QM) or semiclassical propagation, where 
the latter means Vlasov- or Vlasov- Uehling-U hlenbcck 



(VUU) sche mes. Th ese will be reviewed in Sees. III.B.l 



[II. B. 2 and III.B.3. Very violent processes exceed the 



capability of DFT methods and are treated in a purely 
classical manner, either with molecular dynamics (MD) 
or, more simple, with rate equations. We will bricfiy 



sketch both methods in Sees. [II. B. 4 and [II. C 



The large ionic mass usually permits their classical 
propagation by MD. This may be performed simultane- 
ously with the (nonadiabatic) electron cloud or in Born- 
Oppenheimer (BO) approximation, if the electrons follow 
adiabatically the ion field. Light elements (particularly 
H and He) often call for a quantum mechanical treatment 
also for the ions. A full quantum treatment for both, all 
electrons and ions, is extremely demanding and has not 
yet been applied to clusters. However, a QM trea tment 
of He atoms has been widely used for He clusters ( Serra 



et al, 1991; Weisgerber and Rcinhard, 1992 ) and for He 
mater i al in contact with metal clust ers ( [Ancilotto and 



Togi4 |1995| ; [Nakatsukasa et al.[ |00|) 

Figure |^ complements Tab. || in sketching the regimes 
of applicability of theoretical models in the plane of exci- 
tation and particle number. As the decision for a method 
depends on several other aspects (e.g., demand on preci- 
sion, material, time span of simulation), the boundaries 
of the regimes are to be understood as very soft with 
large zones of overlap between the models. Note also the 
two intensity scales on top in Fig. ^, which indicate that 
limitations arc also sensitive to the nature of the system 
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response, i.e., resonant or nonrcsonant. The distirction 
has to be kept in mind when discussing specific systems. 
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FIG. 6 Schematic view of applicability regimes for dilFerent 
approaches in a landscape of system size vs. excitation e lerg; 
The excitation energy can be loosely related to typica 
intensities in the optical range, as indicated by the int 
scales on top for resonant or nonresonant conditions. 
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laser 
5nsity 



Tlic limitations for CI (and other ab-initio met rods) 
are purely a matter of practicability. Time-depe 
local density approximation (TDLDA) is hmited ii 
tern size for practical reasons and in excitation eherg; 
for physical ones, because of the missing dynamica 
relations from electron-electron collisions. The 
limits of VUU are also of purely practical nature 
the lower limits are principle ones, e.g., the neg. 
of shell effects, tunneling, and interference. The 
holds for MD and rate equations. The upper 
energy and/or laser intensity are given by the on^et 
the relativistic regime, where retardation effects 



lim it; 



the Icouplmg begin to severely influence ttie dyn; 



limit results 



For the particle size, a general upper 
the application of the dipole approximation, which 
cally breaks down beyond some ten thousand atoms 
larger systems the field propagation effects (attenu it 
diffraction, reflection) need to be taken into accouit 



B. Effective microscopic theories 



Since a fully ah initio treatment of cluster dyn< 
is hardly feasible, simplifications are necessary by 
inating details of many-body correlations. This 
rally leads to a description in terms of single-p; 
states which is well manageable and still maintain > 
cial quantum features. The eliminated degrees of 
dom are moved to an effective interaction to be 
the re duced description. This leads into the realm 
DFT dPreizler and Gross], |1990D. TDDFT 



namical extension (Gross et al. 
1984), is widely employed in cluster dynamics 
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TABLE III Composition of the basic energy-density functional for electrons, ions and their coupling E = + £"=°"p + 
The ions are described as classical particles with coordinates R/ and momenta P/, / = 1, . . . , Mon- They correspond to the 
nuclear centers and the deeper lying, inert core electrons. The coupling to the electrons is mediated by pseudopotentials Vf^^ 
which are designed to incorporate also the impact of the core electrons on the active electrons. The Vj'^^ counterweight the 
Coulomb singularity of point charges (see Coulomb coupling term) and install effectively a soft inner charge distribution for the 
ion. We show here for simplicity a local pseudopotential which applies throughout all approaches. Nonlocal versions are often 
used in connection with QM electron wavefunctions. The electrons can be treated at various levels of approximation. The QM 
stage employs single-electron wavefunctions (po where a = 1, A^ei- The semiclassical Vlasov description replaces an orbital 
based treatment by a phase-space function /(r, p). In both cases, the Coulomb exchange term and correlations are approximated 
by effective functionals, usually in local density approximation (LDA) and optionally augmented by a self-interaction correction 
(SIC). The fully classical level treats electrons as point particles with specifically tuned effective interaction potentials, e.g., by 
assuming a charge distribution g(r) having a finite width. The total electronic density p(r) is computed differently when going 
from the QM over Vlasov to the MD approaches. Note that the current j(r) is defined analogously to the density. 



2. Time-dependent density-functional theory 

The time-dependent Kohn-Sham (KS) equations cou- 
pled with ionic MD arc derived by variation of the given 
energy (see Tab. EI]) with respect to the single-electron 
wavefunct ions <~p\^ and to the ionic va riables, for details 
see, e.g., ( Reinhard and Suraud , 2003 ). They read 



SE 1 



(9a) 
(9b) 



Since by far most applications employ the local den- 
sity approximation (LDA), the electronic part is coined 
time-dependent LDA (TDLDA). It is coupled to MD for 
the ions, yielding together TDLDA-MD. This treatment 
where electronic and ionic dynamics is propagated simul- 
taneously is compulsory for strong electronic excitations. 

There are many situations where rather slow ionic mo- 
tion dominates and the electron cloud acquires only very 



little excitation energy. For then, one can switch to the 
adiabatic Born-Oppenheimer (BO) picture: 

dtRi = Vp.-Bbo , dtPi = -Vr.-Bbo • (10b) 

It is assumed that the electronic wavefunctions are always 
relaxed into the (electronic) ground state for the given 
ionic configuration and its energy expectation value pro- 
duces a Born-Oppenheimer energy Ebo wh ich d epends 
effectively only on ionic variables, see Eq. (10a). That 
ionic ene rgy -Ebo is then used in a standard ionic MD, 
see Eq. ( |lOb| ). The method allows one to use larger time 
steps because only the slow ionic motion is to be propa- 
gated. On the other hand, full electronic relaxation takes 
many static steps. It depends very much on the particu- 
lar application whether BO-MD is advantageous or not. 

The stationary limit of T DLD A (electronic part) is ob- 
vious - it is given by Eq. (10a). The situation is more 
involved at the side of the ions. A stationary point is de- 
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fined by 9tP/ = and may be reached by simply follow- 
ing the steepest gradient of the potential field. However, 
the ionic energy landscape is swamped by competing lo- 
cal minima. A straightforward gradient path will end up 
in some minimum, but not easily in the lowest one, i.e., 
the ground state. One needs to employ stochastic meth- 
ods, such as simulated annealing and Monte-Carlo sam- 
pling, to explore the high-dimensional landscape of the 
ionic energy surface, for details see ( press et al\ , 1992). 

The most time consuming part in TDLDA-MD, i.e., 
Eqs. (^, is the electron propagation. There are basically 
two different approaches: Basis expansion or coordinate- 
space grid representation, see Tab. |l| and references 
therein. Basis expansions are more efficient in han- 
dling different length scales, as typical for covalent sys- 
tems. Coordinate-space grids, on the other hand, are 
more adapted for the treatment of highly excited systems 
where electron emission plays a crucial role. In the lat- 
ter, absorbing boundary conditions can easily be imple- 
mented to avoid unphysical backscattering for the anal- 
ysis of photoelectron spectra and angular distributions, 
see, e.g., ( palvayrac et al\ , pOOO| ; |Pohl et al\ , p004h| ). A 
very efficient means to find the electronic ground state 
is the accelerated gradient iteration (p31um et al, 1992). 



Time stepping is usually based on a Taylor expansion of 
the time evolution operator. An efficient alternative is 
the time-splitting method which proceeds by interlaced 



kinetic and poten tial evolution ( |Calvayrac et al\ , |2000 



Fcit et al, 1982). The ionic MD u sually empl o ys th e 



1992) 



Velocity- Verlet-algorithm, see, e.g., ( Press et al. 
Ground state configurations are best searched for by 
stochastic methods as mentioned above. 



3. Semiclassical approaches 

As particle number and excitation energy grow, an 
orbital-based treatment of the electronic degrees of free- 
dom becomes practically unfeasible and further approxi- 
mations have to be made. Less demanding are semiclas- 
sical time-dependent density-functional methods, which 
describe the evolution of the one-body electron phase- 
space distribution or the electron density and average lo- 
cal currents. The price for such simplification is the loss 
of the quantized electronic level structure, interference 
effects, and single electron-hole excitations. However, as 
these contributions become less important for larger sys- 
tems with sufficiently narrow energy levels and high ex- 
citations, semiclassical methods provide a powerful tool 
to explore strongly nonlinear laser-cluster dynamics. 

A semiclassical equation of motion for the one-particle 
electron phase-space density /(r, p) as an approximation 
to quantal mean-field dynamics can be found from the 
well-kn own fi —f expansion , see, e . g., (Bcrtsch and Da£ 
Guijtal |1988|; [Pomps et al\, |l997|; [Fennel et al.j |2004|; 



Fennel and K5hn[ , [2008| ; [Plague et al\ , [200C| ). This, to 



lowest order, yields the Vlasov equation 

^/+£-Vr/-Vp/-VrK=ff(r,i) = 0, (11) 

which is widely used in plasma physics. The effec- 
tive electron mean-field interaction potential Vcsir^t) in 
Eq. (^l|) follows from the variation of the potential en- 
ergy Epot = i^coui + E^r. + E'°^P' + E^l, cf. Tab. 0, 



with respect to the local electron density p{r,t), i.e., by 
Vcfi = SEpot/Sp. Ionic motion may be described in the 
same way as for TDLDA-MD, see Eqs. (^). Quantum ef- 
fects, such as exchange and correlation in LDA, are now 
solely contained in the effective potential and the initial 
conditions for the distribution function. The latter can 
be determine d from the s elf-consis t ent T homas-Fermi 
ground state ( Fcrmi[ , 1928| ; Thomas , 1927 ) according to 
/"(r,p) = 2/{27rKye{p^{r) ~p), where 6 the Heaviside 
function, Pp(r) — y^2'm[fi — Vcs(r)] is the local Fermi 
momentum, and fi the chemical potential. The Thomas- 
Fcrmi- Vlasov dynamics resulting from the propagation 
of the initial distribution /°(r,p) according to Eq. ( [Tl| ) 
constitutes the semiclassical counterpart of TDLDA. 

A generic limitation of mean-field approaches, such 
as TDLDA and Vlasov, is the negligence of electron- 
electron collisions. This deficiency may become signifi- 
cant for strong departure from the ground state because 
of considerably weakened Pauli-blocking. In the semi- 
classical formulation, binary collisions can be incorpo- 



Uhlcnbeck (UU) type (Uehling and Uhlcnbeck 




1933), 


see (Bertsch and Das Gupta, 


198?; 


Calvayrac et at, 


2000 


Kohn et al\, |2008|). This results in the Vlasov-Ue 


hling- 



Uhlcnbeck (VUU) equation 

7^/ + - • Vr/ - Vp/ • \/rV,s{r, t) = /uu 
ot m 



with /uu(r, p)= / dfi d'^pi 



(12) 

p-pi|dcr(6', |p-pi|) 
m AD, 

/p'/pi(i - ./p)(i - /pi) - /p./pi(i - ./p')(i - ho 



The collision term embodies a local gain-loss 
balance for elastic electron-electron scattering 
(p>Pi) ^ (p'jPi') determined by the differential 
cross-section d(T(0, jpreiD/dfi, the local phase-space 
densities /p = ./(r,p), and the Pauli blocking factors 
in parenthesis as functions of the relative phase-space 
occupation for paired spins /p = {27Th)^ fp/2. The 
velocity-dependent scattering cross-section can be cal- 
culated for a screened electron-electron potential using 



standard quantum scattering theory ( Domps et al. 
2000| ; [Kohn et al\ [2008| ). Since the collision term in the 



VUU description vanishes in the ground state because 
of the blocking factors, the Vlasov dynamics is recov- 
ered asymptotically in the limit of weak perturbation. 
Commonly, the Vlasov as well as the VUU equation 
are solved by the test particle method only for valence 
electrons, while core electrons are described by ion 



14 



pscud opotcntials, sec, e . g., (Fennel et ai , 2004; [Giglio 
et a{], |2002| ; |Kohn et al\ , ^OQS^ . 
Further simplificati ons can be deduced from hyd rody- 



namic considerations (Ball et 



1973 



Blocli| , |1933| ), i.e., 



by assuming local equilibrium and a slowly varying irro- 
tational velocity field. In this case, the electronic dynam- 
ics can be solely described by the time-dependent elec- 
tron density p(r,<) and a velocity field v(r,t). The cor- 
respondin g equations of motion follow from a variational 
principle (Domps et ai, 1998b| ), leading to a standard 



involved at the same time, which is the typical situa- 
tion in cluster ionization dynamics in strong fields where 
highly charged ions are produced. Hence, very diff'erent 
sets of scale in terms of distances and energies need to be 
resolved numerically, which quickly becomes prohibitive. 

Presently, the single practical solution to microscopi- 
cally resolve ionization dynamics leading to high atomic 
charge states are classical MD techniques. Numer- 
ous groups have developed corresponding methods over 
the years where quasifree electrons and ions are de- 
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netic energy characterizing the local equilibrium an d the 



interaction energy. The continuity equation Eq. (13a) 
and the Euler equation Eq.( |13b[ ) describe the conserva- 
tion of mass and momentum explicitly, while the equa- 
tion of state is implicit in the self-consistent potentials. 
Analogous to Vcff , Vkin results from variation of the now 
density-dependent internal kinetic energy. Within the 
time-dependent Thomas-Fermi (TDFT) approach, the 
internal kinetic energy is described in Thomas-Fermi ap- 
proximation by yjj^(r,0 = |i(37rV(r,i))2/^ TDTF 
represents the most simple semiclassical time-dependent 
density-functional approach. The reduction to the prop- 



Once inncrionized, electrons are explicitly followed ac- 
cording to classical equations of motion under the influ- 
ence of the laser field and their mutual Coulomb interac- 
tion. A striking advantage of the classical treatment is 
the account of the classical microfield and many-particle 
correlations. Nevertheless, there are some difficulties to 
be circumvented. First, the Coulomb interaction has to 
be regularized in order to restore the stability of the 
classical Coulomb system and to avoid classical electron- 
ion recombination below the atomic energy levels. This 
is usually done by smoothi ng the Coulom b inter action, 
e.g., by inserting a cutoff ( Ditmire et al. , 1998) or by 
attrib u ting a n effective wid t h to th e particle (Bclkaccm 



agaiion of four scalar fields tremendously simplifies the et al, 2006b; Fennel et al, 20071 



numerical treatment, which is particulary appealing for 
the study of large systems. For an ap plication to metal 
clusters see, e.g., ( Domps et al , 19981: ). However, as de- 
formations of the local Fermi sphere are neglected (local 
equilibrium), TDTF is not capable to describe thermal 
excitations or highly nonlinear dynamics. 



The second prob- 
lem concerns the computational costs. Standard MD 
simulations scale with the square of the particles num- 
ber due to the direct treatment of the two-body interac- 
tions. For clusters beyond a few thousands of atoms this 
may easily become prohibitive and more elaborate al- 
gorithms such as hierarchical tree codes or electrostatic 
part i cle- in - cell (PIC) methods can be u sed ( Barnes and 



4. Classical molecular dynamics 



Hut, 1986; Pfalzner and Gibbon, 1996). Such methods 



indeed allo w the treatment of large clusters on sufficiently 
long t i mes (llungreuthmayer et al. 2005 Krishnamurthy 



semiclassical, lies in the tact that they arc ot mean-field 
nature and thus neglect the effect of fluctuations, even 
if thermalization due to electron-electron collisions can 
be accounted for approximately in the semiclassical case. 
While mean-field treatments provide a fully acceptable 
approach for moderately perturbed systems, they can- 
not account for the large microfield fiuctuations arising 
from strong-field laser excitation. Exploring these fluctu- 
ations on a microscopic basis requires the construction of 
a statistical ensemble of possible trajectories, which ex 



200S; Baalmann and Rost, 2003, 2005; Saalmann et al. 



A| basic limitation of DFT treatments, quantum or et a^j , p006|; |Kundu and Baue4 | 2006|; pctrov and Davis 



2006). Another option for describing large clusters, even 
at very high laser intensity including rcla tivistic effects , 



are electromagnetic PIC codes, see, e.g., (Fukuda et al. 
2006; Jungreuthmayer et al, 2004| ). 



ceecs standard mean-field capabilities. However, even if 
the approximate description ot strong-field induced clus- 
ter dynamics with the instantaneous ensemble average 
provided by mean- field DFT methods may be sufficient, 
technical difficulties hamper their application to realis- 
tic systems in this case. The problem arises if energetic 
quasifree electrons and strongly bound electrons become 



Inner ionization can be treated in various nonexplicit 
ways. Since deeply bound electrons are associated to 
large energies and short time scales (typically in the at- 
tosecond domain), they are not propagated ex plicitly 
in most cases. An exception can be found in (Belka- 



com et al , 20064 |b|)- In general, however, statistical 



approaches relying on probabilistic estimates of inner 
ionization are used. Common strategies for describing 
atomic field ionization are the consideration of barrier- 
suppression ionization or the application of tunnel ion- 
ization rates, see Sec. [LB. CoUisional ionization may be 



modeled with the semiempirical Lotz cross sections ( Lotz , 



15 



1967). However, this implies that ionization rates, which 
may be altered by many-particle effects in the systems, 
become a crucial ingredient of the dynamics. 



C. Rate equations and the nanoplasma model 

The last step in the hierarchy of approaches from 
the most microscopic to the most macroscopic ones are 
the rate equations models, which describe the system in 
terms of a limited set of averaged global variables. Their 
time evolution is obtained from a few equations account- 
ing for the major couplings, i.e., the interactions with the 
laser field and the internal electronic and ionic processes. 
Such description is based on a continuum picture and 
thus requires the clusters to be sufficiently large. 

The original formulation of a correspon ding model for 



stron g-field cluster dynamics was done by ( Ditmire et al. 



1996) and is known as the nanoplasma model. This 



name reflects the assumption that rapid inner ioniza- 
tion of clusters exposed to intense laser fields creates a 
strongly charged but quasi-homogeneous plasma. The 
typical cluster size domain for which such picture applies 
is the nanometer range, whence the denomination. The 
assumption of a homogeneous plasma requires clusters of 
sizes larger than the Debye length Xd = \/eo kBT/ie^p) 
of the system. Typical density p ^ 10^"^ cm~'^ and tem- 
perature T ^ 1 keV lead to 5 A. 

The basic dynamical degrees of freedom in the 
nan oplasma model are: Nj the number of ions in charge 



stat e the number of "free" (innerionizcd) electrons, 

ii^int the internal energy of the electron cloud, and R the 
radius of the cluster. The global character of these vari- 
ables implies that ions, electrons, and energy are dis- 
tributed homogeneously in a sphere of radius R. The 
evolution of ion numbers Nj follows the rate equation 



dt 



(14) 



where is the ionization rate for ions in charge 

state Nj accounting for tunneling and impact ionization. 
While tunnel ionization dominates early stages of the 
evolution, coUisional ionization takes the lead at later 
times. The electron number N,, evolves according to 



dN,, 



dNi 



dt 



dQ 
dt ' 



(15) 



where Q is the total net charge of the cluster whose 
change is determined by the integrated net flow through 
the cluster surface. The evolution of the cluster radius R 
is determined by the total pressure 



d'^R Pc+Pk 5 



dt^ 



niUii R ' 



(16) 



which is composed of Coulomb pressure pc due to net 
charge and thermal pressure of the hot electron gas 



(treated as an ideal gas of temperature To and internal 
energy E'int = 3/2A^eA;Te). Here, and denote the 
number density and the mass of the ions. 

The internal energy £^int of the electron cloud follows 



dt 



R dt 



^abs = --^5intlm[e(t^las)], 



(17b) 



due to absorption of electromagnetic energy (Pabs), to 
cooling through global expansion {dR/dt term), to ion- 
ization processes (dNj/dt term), and to energy loss 
by electron flow through the cluster surface (Pioss)- 
Here Ip^ are ionization potentials of ions with charge 
state j. The cycle- averaged heating rate Pabs involves 
the volume V and the internal electric field amplitude 
in a dielectric sphere fint = 3fo/(i)/|2 + e((jj)|, where 
£of{t) is the vacuum laser field envelope. The di- 
electric constant e{uj) is usually taken from the Drude 
model e(u>) = 1 — LUp/uj{Lu + iv) with uj^ — n^e^ / (me cq) 
the plasma (or volume plasmon) frequency and the colli- 
sion frequency v for electron-ion scattering. With these 
assumptions, the cycle-averaged heating rate is equiva- 
lent to the expression given in Eq. (|^) in Sec. [I.C and 
exhibits a resonance when the electronic density fulfills 



3 n"'' , where 



mf.f.f^uP' I is termed crit- 



ical density. This condition reflects the Mie plasmon 



resonance of a neutral spherical particle (Kreibig and 



VoUmer, 1995), see the dis cuss ion in Sees. [I.A| and II.C 
The Eqs. (|l4|) through ( |l7b| ) constitute the dynamics 
of the nanoplasma model. In spite of its simplicity, the 
model contains the basic competing processes in the dy- 
namics of the irradiated cluster in a nanoplasma state. It 
is technically simple, but requires several empirical ingre- 
dients, such as, e.g., the various ionization rates. It also 
involves strong simplifications such as a thermal electron 
distribution, an intensity independent heating rate and 
very crude treatment of space-charge effects and electron 
emission. Nevertheless, it was applied to many experi- 
mental results with some successes and its original for- 
mulation was extended in several respects. 

The original model may be questioned at various 
places, e.g., regarding to the assumption of homogeneous 
distributions of all spe cies in the cluster. T his con- 
straint was relaxed in ( [Milchberg et al , ^OOlD by con- 
sidering a radius-dependent distribution. Further, the 
damping effect of the cluster surface is neglected. It 
can be introduced by using a modified collision frequency 
= i^ei + Av / R, which then contains electron- ion colli- 
sions through (vei) and an additional term for surface- 
induced Landau damping Av/R [v - average electron 
velocity). The surface contribution has been shown to 



play an important role for energy absorption (Megi et al. 
200q ). More recently, detailed cross-sections were com- 



puted to include high-order ionization transitions involv- 
ing intermediate excited states for describing the x-ray 
emission from Ar clusters (Micheau et al., 2007] ). An- 
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other important contribution is the lowering of ioniza- 
tion thresholds in the cluster due to plasma screening 
effects (Gets and Krainov. 2006| ), which was shown to 
significantly alter the ion charge d istrib ution as well as 



the heating dynamics (Hilse et al, 2009) 



One should further remind that the nanoplasma model, 
as a statistical continuum picture, may only describe the 
gross features of the interaction of intense lasers with 
clusters. In particular, it cannot access experimental re- 
sults beyond average values. The model may thus fail 
in describing the profiles or far tails of, e.g., ion charge 
state or energy distributions. More detailed insight can, 
for example, be gained from MD simulations. Nonethe- 
less , even in its crudest version the nanoplasma model 



sigh ;s into the time evolution of charging or the explosion 
dynamics for large (nanometer) clusters. 



IV. EXPERIMENTAL METHODS 

With the modern molecular beam machines, the va- 
riety of radiation sources from the infrared to the x-ray 
regime, and the multiply parallel detection and data pro- 
cessing possibilities, challenging and highly sophisticated 
experiments on clusters can be performed. It is possible 
to prepare targets with narrow size distribution or even 
completely size-selected, partially at low or ultralow tem- 



e-g-, 
lani 
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Echt and Recknage^ 1991 




Haberlandl, 


1994 




Mi- 


and lanotta. 


1999; 


Pauly, 


200C; 


Whaley and Millci, 


2001 


). Optical single or many-electron excitation, in 



some cases also being followed by a probing ultrashort 
light pulse, has led to far-reaching insight into fundamen- 
tal processes of the light-matter interaction in clusters. 
In this chapter, rather than covering the vast multitude of 
experimental methods, we review selected current tech- 
niques used for probing dynamics on free clusters. 



A. Generation of cluster beams 

Rare-gas or molecular clusters are produced from an 
adiabatic expansion through a continuously working or 
pulsed nozzle with nozzle diameters ranging from a few 
to 500 ^m, usually restricted by the pumping speed of 
the apparatus. Mixed clusters are generated by a co- 
expansion of a gas mixture or by using a pick-up tech- 
nique with a cross-jet. The cluster size may be varied 
by changing the nozzle temperature or the stagnation 
pressure. Typically, the width AA^ (FWHM) of the size 
distribution roughly equates the average number (N) of 
atoms per cluster. According to semiempirical scaling 
laws (Hagena, 1974, 1981| , 1987) derived from general 
considerations about condensation kinetics, (A^) scales 
with the "condensation parameter" 



r* = fc 



(po /mbar d / /im) 



0.85 



may serve as an acceptable starting point to get first in- et al. ,_ 1996 



where pq is the stagnation pressure, Tq is the nozzle tem- 
perature, d is the effective nozzle diameter. The gas con- 
stants k (in units of [K^-^^^'^mbar~^]) can be calculated 
from the molar enthalpy at zer o tempe r ature and the 
density of the solid according to (Hagena, 1987), ranging 
from k = 185 for Ne, over k = 1646 for Ar and k = 2980 
for Kr, to k = 5554 for Xe. Equation ( [l^ ) holds for 
monoatomic gases; otherwise the exponents of d and Tq 
are different. For conical nozzles, d has to be replaced by 
an equivalent diameter that depends on the half open- 
ing cone angle. The scaling laws developed for rare gases 
have been modified afterwards for metal vapors. 

For experiments at ultra-low temperatures, helium 
drople t pick-up sources prov e to be v ery versatile ([Bartelt 



Goyal et al.j |1992| ; ^Tiggesbaumker and 



Stienkemeiei , 2007). A sketch of a typical setup is shown 
in Fig. 0. He droplets are produced by the supersonic 
expansion of precooled helium gas with a stagnation 
pressure of 20 bar through a 5/im diameter nozzle. By 
choosing the temperature at the orifice (9-16 K), the log- 
normal droplet size distributions can be adjusted in the 
range of (N) = lO'^-lO^ atoms. After passing differential 
pumping stages the beam enters the pickup chamber con- 
taining a gas target or a heated oven, where atoms are 
collected and aggregate to clusters inside the He droplets. 
With this setup it is possible to re c ord cl usters with up to 
150 silver at oms (Radcliffc et al. , 2004) or 2500 magne- 
sium atoms ( Dicderich et al , p005 ), respectively. Down- 



{To/K 



,2.2875 



(18) 



stream another differential pumping stage, laser light or 
an electron beam ionizes the doped droplets. The ben- 
efits of pick-up sources rely on the feasibility to embed 
clusters into a well-controlled environment. In the case 
of He, the embedding medium is superfluid, weakly in- 
teracting, and ultracold with a temperature of about 
0.4 K (Hartmann et al, 1995 ), being an ideal n a noma - 
trix for spectroscopic studies (Whaley and Miller, 2001). 
Similarly, droplets or particles of other elements might 
serve as pick-up medium, e.g., Ar, Kr, or Xe. Subse- 
quent atom agglomeration also can lead to the forma- 
tion of electronically excited species (Eevlev et al., 2000). 
While in particular in the case of helium the nanoma- 
trix is mostly transparent under low laser intensity con- 
ditions, it may become an active part in the interaction 
process under strong laser fields that substantially alters 
the clusters dynamics. Subsequent to plasma formation 
in the embedded cluster, the nanodroplet may be ionized 
as well, giving rise to a core-shell type nanoplasma. 

In these days pure metal clusters are mainly produced 
with laser vaporization or plasma-based methods. In 
both cases the material is vaporized, being partially ion- 
ized, and then undergoes cooling and expansion in a 
rare gas. This can be pulsed, allowing for a hard ex- 
pansion of the seeded clusters into vacuum, or continu- 
ously streaming at lower pressure. In a laser vaporiza- 
tion cluster source a rotating target rod or plate of the 
desired material is mounted close to a piezo or magneti- 
cally driven pulsed gas valve. Usually He pulses with an 
admixture of Ne or Ar at backing pressures of 2-20 bar 
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FIG. 7 Schematics of a He droplet pick-up cluster beam ma- 
chine. Atoms from the vapor in t he pick-up cell can be l oaded 
into the droplets at 0.4 K. After ( Dicdcrich et al , 2005 ) 



serve as seeding gas. Intense ns laser light pulses with 
about 50 to lOOmJ/pulse erode target material by pro- 
ducing a plasma plume, which is flushed by the seeding 
gas through an about 1 mm diameter channel and a noz- 
zle into high vacuum. The close contact with the cold gas 
leads to supersaturation and eflicient aggregation already 
in the source channel. The nozzle - often elongated by 
an extender - can be cone-shaped or merely be a cylin- 
der. In some cases an additional small mixing chamber 
between source body and extender might increase the in- 
tens itv within a desired mass range. Depending on ma- 



teri al and operation conditions, different types of nozzles 
are in use, partially with very long extenders of 10 cm or 
more. There is no optimal photon energy, but the inten- 
sity must be sufficient to induce vaporization or create a 
plasma. However, frequency doubled NdiYAG laser are 
often used, as its green color facilitates the beam ad- 
justment. With laser vaporization sources practically all 
solid materials can be vaporized. As a significant fraction 
(~ 10%) of the clusters is charged, no additional ioniza- 
tion is necessary for studies on mass selected species. 

Several types of plasma-based sources are commonly 
used, the most prominent being the magnetron sputtering 
cluster source, going back to developments in the group 
of Haberland (Haberland et al, 1992| ). The basic erosion 



process is high pressure (1 mbar) magnetron sputtering. 
This versatile tool operates with a few cm in diameter 
plane solid target mounted close to an axial permanent 
magnet, see Fig. ||. In the presence of the seeding gas, 
a high voltage between a ring-shaped electrode and the 
target initiates and drives a discharge, efficiently eroding 
the material and producing a circular well after several 
hours of operation. The mainly charged vapor is cooled 
by the seeding gas and transported through a nozzle. 
Conducting materials can be sputtered by this source, 
whereas ferromagnets may cause difficulties. 

In contrast to the magnetron sputtering source which 
operates with a high voltage discharge, arc cluster ion 
sources make use of high current arcs. Such are known 
as vacuum arcs, self-stabilizing at about 40 Volts and 
40 Amps. The discharge can be sustained in vacuum once 
a spark has initially brought some metal into the vapor 
phase. It is important that the discharge is carried by 
the metal vapor rather than by the seeding gas. In or- 
der to accomplish this, the temporal development of the 
high voltage-driven spare needs special care. Once the 
metallic component in the source rules the conductivity. 




FIG. 8 Plasma plume of an uncovered Haberland-type mag- 
netron sputtering cluster source during operation. The ion 
and electron motion is guided by permanent magnets behind 
the target. Right: picture of a used silver target. 



the discharge voltage switches to a low level so that the 
seeding gas will not directly be ionized. Two variants of 
the arc sources are in use, pulsed ones and continuously 
working ones. The concept of the pulsed arc cluster ion 
source PACIS ( |Cha et al\ , |1992| ; [Siekmann eTaI| , |l99lD is 
very similar to the laser vaporization cluster source, only 
that the laser is replaced by a pulsed high-current arc 
between two electrode rods at about 1 mm separation. 
An offspring of the PACIS uses one rotating electrode, 



then being called "Pulsed Microplasma Source" (Bar- 



borini et al , 1999 ). When operated cont inuously we ob- 



tain t he Arc Cluster Ion S ource ACIS (Kleibert et al 



20071 ; iMcthhng et al[ |200l| ). Here the target is a water- 
cooled hollow cathode, a water-cooled counter electrode 
serves as anode. Magnet coils around the hollow cath- 
ode help to control the arc. Again, the plasma is flushed 
by an inert seeding gas into vacuum, producing a cluster 
beam with a high amount of charged species (about 80%, 
depending on the material). The beams from the ACIS 
can be focussed by aerodynamical lens systems. These 
are sets of orifices and/or confining tubes connected to 
the nozzle. By choosing appropriate dimensions the on- 
axis intensities increase, which goes along with a narrow- 
ing of the particle size distribution (Passig et al., 2006). 
This type of source can generate large metal particles 
from 2 to 15 nm in diameter, an interesting size range for 
future studies of the intense laser-cluster interactions. 

All cluster sources described above are housed inside 
well-pumped vacuum chambers in order to reduce the gas 
load at the point of investigation. Ideally, only the central 
filament of the jet passes a narrow skimmer and enters 
as collimated cluster beam the photoexcitation chamber. 
Further differential pumping can lead to sufficiently low 
pressure for the spectroscopy on isolated species. How- 
ever, many strong-field experiments do not make use of 
single cluster excitation. In particular for rare gas clus- 
ters, the laser is often focussed onto the beam in the high 
pressure zone close to the nozzle. In such cases many 
interacting clusters are simultaneously excited, thus the 
observed signal might originate from a dense cluster en- 
semble rather than from isolated systems. 
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B. Sources for intense radiation 



Within the last 20 years uhrashort-pulse lasers have 
undergone dramatic improvements with respect to pulse 
width, power, and repetition rate. This was first enabled 
by the technique of colliding pulse modelocking (CPM) 



within a ring dye laser (Fork et ai, 1981) and later by 



the invention of the chirped pulse amplification (CPA) 
scheme by Maine et al. (Maine et ai, 1988). Nowadays, 



the broadband fluorescent (690-1050 nm) laser crystal 
Ti:sapphire operating at a central wavelength of 800 nm 
is the working horse in delivering ultrashort and intense 
optical radiation. Laser pulse dm'ations as short as some 
femtoseconds ( Brabec and Kraus: 
attoseconds ( Corkum and Kraus: 



200 



OOOt [Kelleij , g003|) or 
as well as pulse 



powers in the Petawatt regime (Ledingham et ai, 2003) 
are available. To avoid damage of the optical compo- 
nents, the pulses from a modelocked femtosecond laser 
oscillator is first stretched to some ps before amplification 



and then re-compressed in the final step (Maine et al. 



1988). For energy enhancement regenerative amplifiers 



or bow-tie shaped multi-pass configurations are typically 
used. Stretching as well as compression of the pulse is 
achieved by introducing diffractivc elements, e.g., ve^iec- 
tion^trratinfys ([Strickland and Mouroii |l 98."i|'l in the opti- 



cal ; )ath. High energy pulses in other wavelength regions 



can Ibe realized, e.g.. bv amplification of the third har 



monic in a KrF amplifier operating at 248 nm (Bouma 
et a[], |1993 ). Due to the limited bandwidth of the tran- 
sition the pulse duration in this type of laser is limited 
to some hundreds of femtoseconds. With high harmonics 
(HH) generated by focussing intense pulses into atomic 
gases the short wavelength regime becomes access ible 
ope ning up the route t owards attosecond pulses (Pa- 



padogiannis et ai, 199£). Pulse intensities as high as 



1.3 X 10J±W/_ cm have been r eported for the 27th har- 
monic ( Nabekawa et ai , 2005 ). Only recently the vac- 
uum ultraviolet free electron laser (VUV-FEL) FLASH 
at DESY has been setup, currently delivering pulses 
with wave lengths down to 6.5 n m at peak energies up 



to 100 ( lAyvazyan et alj |2006| ) 



In the optical domain single-shot aut ocorrelators or 
more sophisticated setups (Trebino, 2002) are applied for 
pulse characterization. In many experiments only the 
pulse width is varied by detuning t he co mpressor length. 
This introduces a linear chirp (Sec. II.B ) and allows con- 
tinuous variation of the pulse duration between sub-lOOfs 
to many ps. To generate dual-pulses with variable opti- 
cal delay (pump-probe) the initial pulse may be split into 
two replica, e.g., by a Mach-Zehnder setup. Moreover, 
liquid crystal spatial light modulators, acousto optical 
modulators, and deformable mirrors allow one to modify 
the pulse structure at will (Weinei, 2000| ). Besides am- 
plitude and phase, also the polarization can be altered, 
e.g., to drive reactions selectively into a desired channe l 
in coherent control experiments (Brumer and Shapirc , 
1995; Tannor et ai, [1986|). This scheme conn ected to a 



feedback algorithm (Judson and Rabitz, 1992) is capable 



of optimizin g the lase r-matt er coupling, see e.g. (Assion 



et al\ , |1998| ) and Sec. |VII.A 



For pulse focussing, lenses or parabolic mirrors can be 
used. The latter avoids pulse modification due to the 
propagation through optical elements, i.e., pulse broad- 
ening, self-focussing, or phase modulation. The waist 
radius of a Gaussian beam at the focus is wq = 2Xf /n, 
where the /-number relates the size of the unfocussed 
beam diameter D to the focal length of the lens df by 
/ = D/df, and A is the wavelength. Typical spot sizes 
are a few tens of /im. For a qualitative description of 
nonlinear laser-matter interactions the intensity profile 
in the focal region has to be taken into account. For 
a given p eak intensity /p, the inten sity profile /(r, z) is 
given by (Milonni and Eberly, 1988) 



/(r,0) 



1 + z2/z2 



exp 



2r' 



(19) 



where r and z are the axial and transverse distances to 
the focus and zq = ttwq/X specifics the Rayleigh length, 
where the beam radius has increased to ^/^wq ■ The focal 
intensity profile leads to volumetric weighting, which has 
been used to determine in tensity thresholds i n the strong 



et al 



et al. 



2006| ; iGoodworth 



field io nizat i on of atoms (|B r yan et al., ^ , 

2005| ; Hansch et al. , 199E ) and molecules (p3cnis 
2004|). Applied to clusters, this intensity-selective 



scanning method has revealed a dramatic lowering of the 
threshold intensities for p roducing highly charged ion s 
when compared to atoms (Doppner et al., 2007b, 2009). 



C. Particle detection techniques 

Optical excitation of clusters can lead to extensive 
fragmentation. Usually fragment mass spectra are ana- 
lyzed i n terms of stabil i ties, s imilar to nuclear fission pro- 
cesses ( ^chmidt et al , 1992 ). In strong fields, however, 
dedicated techniques are needed to resolve the emission 
spectra of ions and electrons in detail. 



1. Determination of charge state distributions 

The most straightforward method to determine charge 
state distributions of clusters and their fragments is 
ion mass spectrometry. Irrespective of the particular 
method, the mass separation will always be connected 
to the charge-to-mass ratio. In particular time-of-flight 
(TOF) methods with accelerating electrical fields are 
widely used for analyzing charged products after pho- 
toionization. Fig. |9| shows an example of highly charged 
atomic ions emerging from silver clusters embedded in He 
droplets after irradiation with intense fs laser light. The 
TOF spectrum exhibits contributions of He and Ag clus- 
ters with high masses (not shown here). At short flight 
times a situation appears like in Fig. ^. Whereas the 
background peaks are signatures of the He droplet frag- 
ments, the highlighted series can uniquely be assigned to 



19 



atomic ions in high charge states from the Coulomb ex- 
plosion of Agjv- As a matter of fact, the Ag ions carry 
high recoil energies due to the violent expansion. There- 
fore TOF methods that use an acceleration of the ionic 
ensemble by electric fields in the few kV range loose part 
of their resolution and transmission. Consequently, the 
TOF spectra only prove the occurrence of the ions but 
do usually not image the real charge state distribution. 



2. Acquisition of ion recoil energy spectra 

A simple and versatile tool to investigate ion recoil en- 
ergies is the acceleration-free TOF spectroscopy. Two 
preconditions have to be met in order to allow a unique 
interpretation of the results: First, there has to be a 
defined source point for the ion emission. Second, the 
nature (mass) of the ions must be known, which often 
is a point difficult to achieve. However, the excitation of 
single-element clusters with sufficiently strong laser fields 
leads to complete fragmentation into atomic ions with 
known mass. In this case, the kinetic energy is deter- 
mined by TOF measurements through a field free drift 
tube of about 0.5 m, without initial electric field. For 
reducing noise caused by secondary electrons and, more- 
over, to restrict the ion detection to the Rayleigh region 
of the laser focus, an adjustable narrow slit confines the 
ion trajectories. Resulting TOF spectra can then directly 



be c mverted into kinetic energy spectra, sec, e.g.. Fig. l|d. 

The field-free ion TOF yields recoil energies irrespec- 
tive of the ion charge states. For a detailed analysis 
it is necessary to resolve charge state dependent recoil 
energies. To this end two methods have successfully 
been applied, both of which simultaneously measure the 
ion charge state and energy. The first one uses mag- 
netic deflection time-of-f light (MD-TOF) mass spectrom- 
etry ( Lezius et ai , 199§| ). This technique bases on TOF 



measurements at different positions behind a magnetic 
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FIG. 9 Charge state spectrum from a time-of-flight analysis 
of Agjv in He droplets with (N) = 40, exposed to 400 fs laser 
pulses at 4 X 10"W/cm^ and 800 nm. The resulting Ag''+ 
signals from the Coulomb explosion are highlighted. Ions with 
up to (jf = 11 are been detected. The occurrence of He^''" 
stems from charge tr ansfer with the Ag io ns at the chosen 
laser intensity. From (Doppner et ai, 2005). 



field. With the MD-TOF, highly energetic (up to 1 MeV), 
multiply charged ions could be recorded. 

The second method is of static nature and bases on 



a principle first applied by Thomson ( Thomson , 1907 ). 
Fig. |l^ sketches the Thomson analyzer for the simulta- 
neous measurement of energy and charge of ions expelled 
from an exploding cluster. It consists of parallel electric 
and magnetic fields, followed by a field free drift zone 
in connection with a position-sensitive detector. The ex- 
perimentally obtained raw data reflect momentum and 
energy per charge and have to be transformed to en- 
ergy vs. charge spectra. For Agw the charge state re- 
solved ion energy distribution is rather narrow and the 
maxi mum energy grows alm ost linearly with ionization 



stage ( Ddppncr et al , 2003 ) 



3. Energy and angular resolved electron detection 

The experimental challenge in photoelectron spec- 
troscopy results from the notoriously low densities in 
mass-selected charged cluster beams. To cope with this, 
time-of-flight electron spectroscopy has been developed 
with a magnetic field gradient. When the clusters are 
ionized at a certain spot within an electron magnetic 
bottle spectrometer the complete photoelectron spectrum 
can be recorded by time-of-fli ght measurem e nts w i th up 



to 100% de t ection efficiency ( Arnold et al j_ 1991 Gan- 
tefor et ai , 1988 ; Kruit and Read , [1983 ; Taylor et al. 



1992). Whereas this method turned out to be extremely 



fruitful to reveal the electronic level structure of many 
mass selected cluster anions, the magnetic fields involved 
hamper the retrieval of satisfying angular information. 
In the case of a neutral cluster beam, the target den- 
sity can be sufficiently high in order to get a spectrum 
even without the magnetic field. Electron emission and 
drift occur within a field-free tube, equipped with a time- 
resolving detector. By rotating the polarization direc- 
tion of the laser, angular resolved photoelectron spectra 
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FIG. 10 Sketch of the Thomson analyzer. Ions enter a region 
of parallel electric and magnetic fields trough a tiny hole. The 
resulting defiection gives characteristic parabolas from which 
the charge state selective recoil energy can be deduced. A 
multi channel plate detect or with an imag i ng sy stem serves 
to record the data. From (Doppner et al, 2003), with kind 
permission of The European Physical Journal (EPJ). 
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are obtained. An increasing length of the drift tube in- 
creases the energy resolution on the expense of signal 
intensity. Acceptable results can be achieved with mag- 
netically well-shielded tube of about 0.5 m length. 

In contrast to the electron TOF method, where kinetic 
energy release information is contained in the electron 
drift times, imaging techniques extract energy and angu- 
lar distributions from spatially resolving detection. The 
striking advantage of this method is that the full emission 
characteristics can be reconstructed from the 2-D image 
by means of an Abel inversion. The energy resolution 
is limited by the quality of the 2-D detector (Heck and 
1995| ). An improvement of the 2-D imaging 



Cha iidler 



technique has been obtained by introducing a lens op- 
tics which maps all particles with the same initi al vcloc- 
ity vector o n to th e same point on the detector (Eppink 
and Parker , 1997 ). So far, this technique has mainly 
been used to record low-energy electron spectra. With 
modified electrode configuration energetic electrons from 
clust ers driven to Coulomb e xplosion are accessible as 
well (Skruszewicz et al, 2009). 



V. SINGLE- AND MULTIPHOTON PROCESSES IN 
CLUSTERS 

The previous sections have provided basic tools for the 
description and analysis of laser-induced cluster dynam- 
ics. In the following presentation of specif ic exa mples we 
begin with single-photon process es in Sec. V.A and move 
on to multiphoton effects in Sec. V.E| . In both cases clear 
signatures of the photon energy persist. Single-photon 
excitations are typically investigated by photoelectron 
spectroscopy (PES), which is usually interpreted as a 
static image of the density of states and so indirectly of 
the underlying geometry. When carried out with angular 
resolution, PES reveals structural details of the electronic 
orbitals being excited. However, even single-photon pho- 
toemission goes beyond a mapping of system properties 
in a static and direct way, as it reflects a dynamical pro- 
cess. Pump-probe studies, as a time-resolved version of 
PES, give access to ultrafast structural dynamics and en- 
ergy redistribution pathways. Additional reaction chan- 
nels emerge with the abso rptio n of multiple photons, as 
will be the subject of Sec. V.B[ Besides above-threshold 
ionization, as a prime example for multiphoton signa- 
tures, thermalization and its effect on electron spectra 
will be discussed. Another issue are plasmons. which 
often govern the response of metal clusters and become 
broadened by nonlinear contributions at higher intensity. 
However, they remain a dominant doorway process up to 
the strong-field domain, which is subject of Sec. Vl. 
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FIG. 11 PES spectrum of Na^j^: experimental result (black 
curve) from nanosecond laser excitation with hcj = 4.02 eV 
at T ~ 100 K and theoretical DOS calculated by DFT us- 
ing different ground state structures (as shown). From the 
matching of the spectra the left structu re is favored while t he 
right ones show less agreement. After (Kostko et al, 2007) 



A. Single-photon electron emission 

1. Probing the density of states 

For studying single-electron excitations by photoemis- 
sion it is often u seful to ass u me, m otivated by Koop- 
mann's theorem (Weissbluth, 197S), that the essential 
structures of the electron and ionic systems do not change 
significantly upon electron emission. The photoelectron 
energy spectrum thus basically images the density of 
states (DOS). Based on this assumption, PES has be- 
come a powerful tool to explore the electronic structure 
of mass-selected clusters. Figure ^ displays an exam- 
ple from Na^]^. The measured data (black curve) ex- 
hibits pronounced peaks at binding energies between 1.8 
and 3.5 cV. Such electronic fingerprints reveal details of 
the quantum confinement and change dramatically with 
cluster size or structure. With DFT calculations it has 
become possible to obtain theoretical DOS for compar- 
ison with experimental PES spectra. Fig. |ri| displays 
an attempt to identify the cluster ground state geometry 
out of theoretically suggested candidates by matching the 
DOS. A vast amount of photoelectron spectra on differ- 
ent systems has been accumulated since first successful 
experiments (Cheshnovsky et al.[ |1990|; iGantefor et al 



1987; McHugh et al, 1989; Pettiettc et al, 1988). Dur- 



experiments ( unesnnovsKy et at. ^ Lyyu uanteior et at. 
19961 ; [Gantcfor et aL| |1984;|Ho et al\, |1990|; [Leopold et al 



ing the course of time, developments in cluster produc- 
tion and electron detection have made it possible to cover 
large size ranges at high energy resolution. For instance. 



in ( Wriggc et al , 2002 ) PES spectra of Na^ for N=31-500 



show peaks that can be assigned to the electronic shell 
structure. For small systems a higher level of theoretical 
understanding can be obtained from ab-ini t io qu antum 



chemical methods ( Bonacic-Koutccky et al, 1991) 



To date, most PES studies rely on low-energy pho- 
ton excitations, i.e., valence-band PES. Inner shell pho- 
toionization. i.e., core-level PES, has been demonstrated 



as well (Ebcrhardt et al , 199C; Siekmann et 



1993 



Werthcim, 1989). These studies, however, dealt with de- 
posited clusters excited with high photon energy lamps 
or synchrotron radiation. Common results are shifts of 
core levels with cluster size. Due to the surface contact a 
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thorough understanding remains difficult since core-hole 
screening, chemical shifts, electronic relaxation or charge 
transfer dynamics contribute to the spectra. 

With third-generation synchrotron sources, experi- 
ments on free neutral (not mass-selected) clusters became 
possible. One issue of such studies is the absorption site 



as a probe of the local enviro nment (Hatsui et ai, 2005 



von Pietrowski et ai, 2006) 



measured line profiles ( rchaplyguine et al 



In rare gas clusters the 
^004|) show 



well-separated features that can be attributed to the ion- 



ization of_surfe££_aj2d_voI^^ atoms, respectively (Amai 
et a[], ^005 ; Bergersen et al., 2006). Such analyses can 



also provide an indirect size measurement, as has recently 
been shown for neutral nanometer clusters of various 



metal s, i.e.. 



Na ( Peredkov et alj |2007a|), Pb (Peredkov 



et a[], |2007b ), Cu, and Ag ( Tchaplyguinc et al. 



2007) 



Latest progress in core-level PES has been achieved 
at the free electron laser FLASH which delivers intense 
pulses with up to 200 eV photon energy. The energy 
range and high brilliance open new possibilities to in- 
terrogate both the complete valence regions as well as 
shallow core levels of numerous systems. For example, 
PES on free mass separated Pb^ revealed a pronounced 
A^-dependent shift of the bd core level ( penz et~al. 



which is in accordance with the metallic droplet picture 
for large N . However, strong deviations starting below 
A'^ < 20 indicate a transition from metallic to nonmetallic 
bonding due to less efficient core-hole screening. 

A solid theoretical understanding of the photoioniza- 
tion process requires the complete toolbox of computa- 
tional many-particle physics. One example where DFT 
calculations for Na^ are compared to experimental PES 
was shown above in Fig |l^. In the same spirit, Si^ for 
A^ = 20 — 26 have been investigated theoretically in (Gu- 



lian ov et qL|, |2GG5| ) and compared to data from (Hoff- 



marn et 7 |2001^ In both cases, not all peaks could be 



fully reproduced by theory, especially for deeply bound 
electronic states. Nevertheless, from comparison of the 
calculated DOS with the experiment the ground state ge- 
ometry can be identified and discriminated against com- 
peting isomers in many cases. Remaining discrepancies 
reflect that static DFT calculations based on the Kohn- 
Sham eigenvalues are insufficient to fully describe the 
photoemission. It is well-known that the interpretation 
of eigenvalues as single-particle energies requires atten- 

p006| ) 



tion (Kiimmel and Kronik, 2008; Mundt et al. 



This concerns the meaning of single-particle eigenvalues 
itself as well as dynamical aspects, as Koopmann's theo- 
rem does not hold in a strict way. In other words, photoe- 
mission is a highly correlated process. The photoelectron 
interacts with the residual system during its removal and 
may substantially modify the level structure. The effect 
becomes important with low energy electrons and dra- 
matic in the zero electron kinetic energy measurements. 

The question whether PES reflects parent or daugh- 
ter cluster DOS or a dynamical mixture of both has 
been tackled in the case of sodium cluster anions, see 
Fig. n% The comparison between the experimental spec- 
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FIG. 12 Comparison of measu red PES spectra for Naf (lower 



panel) (Moseler et al, 2003) and two different theoretical 
predictions. The upper panel shows the single-electron lev- 
els from a (static) Kohn-Sham calculation applying ADSIC. 
The middle panel presents the theoretical result deduced from 
the excitation spectrum of neutral Nar, the final product 
after pho toemission. The excitatio ns were computed with 



TDLDA ( [Mundt and Kiimmcj , |2007] ) 



trum ( Moseler et al. , 2003 ) and the Kohn-Sham eigenval- 
ues of the {parent) cluster anion calculated with average- 
density self-interaction correction (ADSIC) is clearly not 
satisfying (Legrand et al., 2002). A way to circumvent 



the use of the Kohn-Sham eigenenergies is to perform 
a time-dependent DFT calculation of the re s ponse to a 
small pertubation. In (Mundt and Kiimmel, ^007 ), the 
energies of excited states of the neutralized daughter clus- 
ter are extracted from the time evolution of the dipole 
and quadrupole moments and are related to the photo- 
electron kinetic energies by energy conservation (middle 
panel). While some discrepancies still remain there is a 
clear improvement over mere static considerations which 
points out the key role of final state interactions. 



2. Angular distributions 

Besides pure energy spectra, which reflect the elec- 
tronic level structure, photoemission may also reveal de- 
tails of the involved orbitals and thermalization phenom- 
ena. To this end the emission has to be analyzed with 
angular resolution, a subject that still is in its early stage. 
The directionality of the photoelectron angular distribu- 
tion (PAD) can be quantified by a Legendre expansion: 



do- _ CTtot 

d^ ~ 47r 



\l 



(20) 



where 6 denotes the emission angle with respect to the 
laser polarization axis. The anisotropy parameter (32 
ranges from — 1 (emission perpendicular to polarization) , 
over (isotropic emission), to 2 (emission parallel to po- 
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larization) and depends on the orbital symmetry of the 
initial and final states and the electron kinetic energy. 
Only a few PAD experiments have been performed 



on clusters so far 
iV = 4- 
and Hg^, 



Among t hem arc results on W 



2001 



11 ( Baguenard et al. 

20 (|Verlet et al. 



N ^ 3 



P inare et al.[ |1999D 
20041 ). The corre 



sponding (32 as a function of N arc presented in Fig. 
For WjY, (^2 is changing from a more directed behavior 



with small clusters (/32 
sion {P2 0) when N - 



1) to nearly isotropic emis- 
11. From the rapidly reached 



isotropic behavior it was concluded that larger show 
an indirect emission process, where electron-electron col- 
lisions lead to a loss of coherence. This is in line wit h the 
tendency of to undergo thermionic emission (Leis- 
ner =t al , 1991 ). Figure |l^ further shows results on Hg^ 
with strongly size-dependent asymmetries. Although the 
physical origin of these /32- fluctuations could not be clar- 
ified yet, the data illustrate the high system-sensitivity 
of angular resolved photocmission. 

A clear dependence of the PAD on the electronic 
level being excited was demonstrated with medium-sized 
Na^ (Bartels et al, 2009| ). Exemplarily, Fig. ^ com- 
pares a standard PES spectrum of Na^g (top) with the 
corresponding angular resolved result (bottom panel). 

The peaks in the top panel can easily be attributed 
to emission fr o m th e 2p, Ig and 2d shells, see ( Brackj , 
1993; |dc Heei , 1993) for details on the shell nomencla- 
ture. For the given photon energy, the comparison with 
the PAD shows that the 2p and 2c? electrons are emitted 
parallel to the laser polarization, while the Ig emission is 
aligned perpendicularly. These results demonstrate that 
the cluster valence electrons preserve their angular mo- 
menta after excitation. Further, similar angular distri- 
butions for electronic sublevels of a particular shell (e.g. 
Ig) show that the ionic background does not destroy the 
free angular momentum eigenstate character within an 
electronic shell. Thus, the results justify a single-particle 
picture of almost free delocalized electrons for describing 
the PAD from simple metal clusters. 
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FIG. 13 Anisotropy parameter P2 extracted from photoelec- 
tron angular distributions as a function of cluste r size N : Wj^ 



(squa res) exposed to 4.025 eV laser light, from (Pinare et al 



199!:), with kind permission from Springer ScienceJBusmess 



Medi a; Hg^y (circles) irradiated at 3.15 eV, from ( Verlet et al 
2004), with permission from American Institute of Physics. 
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FIG. 14 Photoemission from Na^g obtained with 2.48 eV pho- 
ton energy. The top panel displays angle-integrated spectra 
while the bottom panel shows the corresponding angular re- 
solved results. The emission a ngle 9 is defined with respect to 



the laser polarization. From (Bartels et al 
with permission from AAAS. 
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FIG. 15 Photoelectron angular distribution of Najj" irradiated 
with three differ ent laser energies as indicated, calculated in 
TDDFT. After (|PohI et all |2004b[). 



Besides a state-sensitivity, PAD spectra are also depen- 
dent on the photon energy. A theoretical study of this 



2004b) and is illus- 



effect was reported in (Pohl et al 

trated in Fig. |l^. The figure shows angular distributions 
for three excitations close to and far above the ionization 
threshold by considering an aligned cluster. The pat- 
terns depend on hujii^s and reflect that the nodal structure 
of the outgoing wave (pk.(v) changes with momentum k. 
Note that the latter is asymptotically related to the ex- 
citation energy by |k| = y^2'm,{hu!\iis — Eip). Systematic 
scanning of ?iwias modulates the zeroes and maxima of ip^ 
and thus in principle allows one to systematically probe 
the orbitals of cluster electrons and, in turn, the back- 
ground ionic field. Therefore, PAD is a promising method 
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FIG. 16 Time-resolved photoelectron spectra ol Agj~9 ob- 
tained with laser pulses of about 3 x 10^ W/cm^, TiWpump ~ 
1.55 eV, and hcOprohc = 3.1 eV. The indicated times are the 
delays between the pump and probe pulses. The dashed arrow 
emphasizes the temp oral development of the initially popu- 
lated If level. After ([Niemietz et all pOO^. 



to gain insight into structural cluster properties. Further, 
as we shall sec below, photoelectron spectroscopy is an 
excellent tool for analyzing dynamical processes. 



3. Time-resolved analysis 

Excited states populated by cluster photoactivation 
can decay in different manners, i.e., by emission of ra- 
diation, by internal conversion with energy transfer to 
the ionic degrees of freedom, or by Auger processes. The 
real time dynamics of such processes can in principle be 
explored by tracing the occupation and spectral positions 
of electronic states within time-resolved PES (TRPES). 
For example, structural changes of the ionic frame might 
open fast radiationless decay channels due to transient 
crossings of the potential energy curves of the excited and 
the ground state (conical intersection). The depletion of 
the excited level, when explored in TRPES experiments 
with fs pump-probe techniques, thus offer insight into in- 
ternal energy conversion processes and ionic relaxation 
time scales. Below we discuss two recent examples. 

TRPES on mass-selected Ag^ with N=3-21 has 
been investigated in a two-color pump-probe experi- 
ment (Niemietz et ai, p007 ), using a 1.55eV pump pho- 
ton safely below the vertical detachment energy (VDE). 
Fig. |l^ depicts a series of spectra from Agj~g obtained for 
time delays up to 5ps. The peak at 1.82 eV observed in 
the top panel is consistent with the value of the VDE 



from the If level previously excited by the pump pulse. 
The feature then shifts with time to lower kinetic ener- 
gies (higher binding energy) and stabilizes between 1.1 
and 1.5 eV after about 1.3 ps (see arrow). Subsequently 
this peak loses intensity and vanishes after about 5 ps. 
This evolution can be explained by a continuous Jahn- 
Teller deformation of the excited cluster which eventually 
opens a nonradiative transition channel to the ground 
state potential. For Agj"g, the decay time is estimated to 
be 630 fs, i.e., much shorter than for radiating transitions. 

Another example on gold cluster anions revealed an ex- 
tremely strong size dependence of the excited state life- 
time. More specifically, Au ^^ shows an exceptio nal long 
lifetime of more than 90ns (Walter et al, 2007). Corre- 
sponding DFT and linear-response TDDFT calculations 
predict decay times of 730 ns. In contrast to that, Au^ 
and Au^ show clear indications for fast internal conver- 
sion on the ps time scale driven by ionic motion. As 
displayed in Fig. p7|a, the excited state peak in the ex- 
perimental data of Auy is initially observed around a 
binding energy of 2eV. This value just reflects the cn- 
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FIG. 17 Time-resolved photoelectron spectra of Auy : (a) Ex- 
perimental data obtained with 40 fs pulses at intensities be- 
low 10" W/cm^ with hij = 1.56/3.12 eV for pump/probe; (b) 
Simulated PES obtained from TDDFT coupled to ensemble- 
MD "on the fly" . Note that the experimental data is plotted 
with a logarithmic del av axis while it is lin ear for the calcu- 



lation results. From (Btanzel et ai, 2007), with permission 



from American Institute of Physics. 

ergy difference between the excited anion and the neu- 
tral ground state. As time evolves, the peak first slightly 
shifts to higher binding energies, then loses intensity, and 
finally becomes washed out after about 1 ps. At the same 
time a new feature appears between 2.4 and 2.8 eV which 
has reached high signal intensity around 3 ps. The time 
evolution of the experimental peak maximum yields an 
exponential decay with a time constant of 1.8 ps. Results 
from a corresponding linear-response TDDFT calculation 
based on a propagation of an ensemble of classical tra- 
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jectories "on the fly" (stanzcl et al, 2007) are depicted 
in Fig. |l7|b. Note the hnear scale used here, whereas 
the experimental data is plotted on a logarithmic scale. 
The calculated population dynamics show a decay time 
of 1.9 ps, which is in good agreement with the experi- 
mental result. A closer analysis of the calculation results 
reveals the following process: The excited anion relaxes 
towards a crossing with the anionic ground state within 
about 340 fs, see the rapid lowering of the initial peak at 
2 eV in Fig. pT^j. At the crossing the excited state begins 
to populate the vibrationally excited ground state. This 
mixing is expressed in the new feature emerging around 
3.0 eV. The rapid bleaching of the 2.3 eV feature beyond 
3ps and the stabilization of the 3.0 eV peak reflect inter- 
nal energy conversion and further indicate a melting of 



the cluster (Stanzel et al, 2007) 



An even more elaborate way for investigating cluster 
dynamics via photoemission is offered by recording time- 
resolved photoelectron angular distributions (TRPAD). 
The applicability has been demonstrated with molecules, 
see (Suzuki, 200(: ). Femtosecond TRPAD studies on clus- 
ters have been performed by the Neumark group using 
imaging techniques to investigate the relaxation dynam- 



ics of small Hg^, see (Bragg et ai, 2005) 



B. Multiphoton signatures 

At intermediate laser intensity, single-photon processes 
as discussed above begin to be accompanied by multi- 
photon effects. The current section is devoted to this 
transition region where linear and nonlinear excitations 
simultaneously occur. We will focus on processes like 
above-threshold ionization, resonance broadening, and 
the onset of electronic thcrmalization. In all cases the 
laser intensity remains sufficiently low to resolve the in- 
fluence of the chosen photon energy. 



1. Competition of linear and nonlinear excitation 

Both single- and multiphoton processes can be demon- 
strated with the exam ple of photoe l ectro n spectroscopy 
on Hgj"4, see Fig. |l| ( [Verlct "alt |00^. These spec- 
tra are extracted from PAD measurements, see the polar 
plots in panels a, b and c and the corresponding angle- 
integrated results at the bottom. Mechanisms leading 
to specific peaks are schematically indicated in the top 
panels. Process A depicts the direct emission of the ex- 
tra electron in the 6p level with one high-energy photon 
or with two low-energy photons. No excitation energy is 
transferred to other decay channels and the peak labeled 
A reflects the energy of the photons and the VDE of the 
6p electron. Other pathways leading to electron emis- 
sion arc feasible if the photon energy exceeds the s — p 
band gap: single-photon detachment of an electron in the 
6s band (case E); two-photon intcrband excitation of an 
electron from the s band via the p band (processes C and 



D); or single-photon Auger processes, labeled B. The lat- 
ter excitation scheme is observed in the spectra as broad 
features with an onset consistent with the VDE of the 
s and the p band in Hgj"^. This example thus demon- 
strates the complex pathways in a seemingly simple case 
of laser-exposed clusters. 

The scenarios exemplified in Fig. |l^ call for a more 
detailed theoretical analysis. Some studies in the sim- 
pler case of Na clusters were performed in a scries of pa- 
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persin the framework of TDLDA-MD ( |Pohl et 
2004a). For our purpose, we concentrate on the impact 



of intermediate states in multiphoton induced photoe- 
mission (Pohl et ai, 2001). In this laser intensity regime, 
the electron single particle energy Sq is a priori deduced 
from the recorded electron kinetic energy by writing 



las ; 



(21) 



where v is the number of photons involved. Chasing such 
a multiphoton process is in principle straightforward as 
the signal shifts according to the change in the photon 
energy. The picture becomes more complicated if an in- 
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FIG. 18 Top: Sketch of possible electron emission pathways 
from negatively charged mercury clusters. Note that before 
irradiation the additional electron is located above the band 
gap (BG) of the corresponding neutral system. Bottom: An- 
gular resolved and integrated photoelectron spectra of Hg]^4 
at different photon energies, (a) and (d) with ficJias = 1.57 
eV, (b) and (e) with ?iwias = 3.15 eV, and (c) and (f) with 
fttjjias = 4.58 eV. The labels cor respond to the pro cesses 
sketched in the top panel. From (Verlet et at, 2004), with 
permission from American Institute of Physics. 
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termediate state (say, of energy ei) can be populated 
by a one or multiphoton process from an initial state 
So- Then the original direct i^-order process competes 
with a sequential v — 1 (or v — 2, . . . ,) process from state 
El. Since sequential processes are less shifted than direct 
ones, discrimination is again possible by slightly changing 
the photon energy. A somewhat similar situation occurs 
when the photon energy is close to the plasmon excita- 
tion (Pohl et al, 2001). In this case the spectra exhibit 
contributions pinned to the resonance energy. 

The plasmon plays a major role in the optical excita- 
tion of simple metal clusters, see, e.g.. Fig. 0b for ex- 
perimental examples. Even at higher laser intensities 
where details of level spectroscopy go lost the plasmon 
remains quite robust. This is exemplarily illustrated in 
Fig. |l^ for NaJi^ and irradiation with about 10^^ W/cm^ 
(top panel) . The plasmon peak centers around 3 eV 
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FIG. 19 Top: Comparison of the total valence electron 
emission from NaJi calculated with TDLDA, Vlasov, and 
VUU calculations as function of photon energy for ex- 
citation with a 20 f s Gaussian l aser pulses of intensity 



7 = 6xlO"W/cm^ (Giglio et al, 2003). Bottom: TDLDA 



results of the valence electron emission from Nag over pho- 
ton energy for excitation with 10 fs pulses for three different 
intensities (as indicated). From (Ullrich et al., 1997). 



and a comparison between various calculations is per- 
formed. The high number of emitted electrons calls 
for calculations going beyond the mere mean-field and 
thus requires the inclusion of dynamical correlations and 
electron-electron collisions, as is done, e.g., with VUU. 
However, the strong impact of the plasmon is obvious 
in all models. As it provides a resonant coupling chan- 
nel, enhanced energy absorption and increased electron 



emission are observed (Calvayrac et al., 200C). It should 
be noted that TDLDA-MD and Vlasov-LDA-MD match 
almost perfectly, but for details in the tail of the distribu- 
tions (e.g. at 3.5 eV). The peak height for VUU is lower 
than for the pure mean-field approaches due to damp- 
ing of the resonance by electron-electron collisions. This 
damping (or resonanc e broad ening) leads to a higher ab- 
sorption (Kohn et al, 2008) and slightly enhances ion- 



ization for off-resonant excitation. 

The nonlinear nature of the response can be probed 
directly by varying the laser intensity, see the bottom 
panel of Fig. |9]. Although the calculations on Na^ have 
been restricted to the TDLDA level (no dynamical corre- 
lations), they qualitatively reflect major trends. In gen- 
eral, the yield is maximal for photon energies in the vicin- 
ity of the plasmon. The peak height strongly increases 
with laser intensity, but docs not show a linear scaling - 
a clear sign of nonlinearity. In addition, the shape of the 
spectrum is substantially altered at higher laser intensity. 
The peak width increases dramatically from about 0.3 eV 
at 10" W/cm2, 0.6eV at lO^^ w/cm^ to almost 2eV 
at 10^'^ W/cm^. This indicates clearly a transition from 
the frequency dominated domain to a dynamical regime 
where spectra become less sensitive to the frequency of 
the laser to the benefit of the field intensity. 



2. Above-threshold ionization and thermalization 

At sufficiently high photon density, two and more pho- 
tons can cooperate almost simultaneously in the exci- 
tation of a single electron. One of the consequences 
is direct electron emission, even far beyond the ioniza- 
tion threshold, although the photon energy stays below 
the IP. Fig. gO| shows computed above-threshold ioniza- 
tion (ATI) spectra for Na^, where the photon energy of 
2.7 eV is to be compared with Eip of 7.5 eV [Ip state) 
and the binding energy of the Is state of 8.8 eV. The 
spectrum at the lowest intensity /q shows distinct peaks 
which can be associated with emission from these states 
and a well-defined number of photons v, as defined in 
Eq. (21). The intensity /q is already intermediate, as a 
sufficient photon density is required to drive multipho- 
ton excitation. On the other hand, the emission from a 
I'-photon process evolves as I'^ . This produces a steep in- 
crease with intensity such that there remains only a small 
intensity window before the signal becomes blurred. This 
is demonstrated in Fig. |2^ with the two higher intensi- 
ties increased by factors of 1.5 and 3.0, respectively. Note 
that the spectrum is almost structureless for the highest 
intensity, which can be explained in the following way: 
The ongoing ionization increases the bonding and down- 
shifts the single-particle levels. Moreover, the spectrum 
collects electrons from all stages such that the level mo- 
tion first induces a broadening of the peaks (7=1.5 /q) 
and finally a complete blurring. It should be noted that 
the fully smoothed distribution shows an exponential de- 
crease (after removing the V^'kin phase-space element). 
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In fact, as a result of the /"-law, the exponential trend 
already appears for lower intensity when connecting in- 



2004a) 



dividual peaks, for details see (Pohl et al 

The example in Fig. |2^ with fixed pulse length deals 
with the intensity effect on ATI from clusters. How- 
ever, the pulse duration plays a central role as well be- 
cause competing perturbation by electron-electron colli- 
sions and ionic motion come into play with increasing in- 
teraction time. This aspect is worked out in Fig. ^ with 
an experimental ATI result on Cgo- The two middle pan- 
els corroborate the previous observation that increasing 
the laser intensity (step from right middle to left middle 
panel) smears out the multiphoton peak structure. Going 
through the figure from top to bottom, i.e., along increas- 
ing pulse duration, one also obtains a disappearance of 
the detailed pattern, but this time due to an increase in 
the pulse width. A quick glance back on Fig. ^ helps to 
an interpretation. After a time span of the order of the 
electron-electron collision time (rge), the photon energy 
is distributed over the whole electronic system. The thus 
thermalized cloud evaporates one or more electrons at 
later times. The corresponding thermal emission spec- 
trum is a smooth exponential (times phase-space factor), 
as indicated by the dashed curves in Fig. |2^. The inter- 
mediate pulse width of 110 fs obviously excites a transi- 
tional stage where direct and thermal emission compete. 
The longer time of 500 fs is safely in the thermal regime 
and even lowering of the intensity (lower right panel) does 
not revive any detailed ATI structures. 

Since pronounced ATI signatures can be observed with 
complex systems like Ceo, it is a logical next step to 
ask for the angular electron distributions as in the case 
of single-photon excitations discussed before in Fig. |l^. 
In a corresponding experiment on Cgo the energy and 
angular resolved spectra where measured for excitation 




FIG. 20 Photoelectron spectra from Na,^ exposed to 100 fs 
laser pulses with fiuj = 2.7 eV at d ifferent intensities (indi- 



cated) as calculated with TDLDA ( |PohI et al\ , |2000|). Ver- 
tical dotted lines show energies corresponding to Eip plus a 
certain number of photons as indicated. Bound states from 
which electrons originate are indicated in a few cases. 



with 60 fs laser pulses, see Fig. Integration of the 
signal over the emission angle basically reproduces the 
sequence of ATI peaks, as seen for the short pulses in 
Fig. ^ The PAD exhibits as an additional information 
that the ATI peaks are well collimatcd in the direction 
of the laser polarization, an effect which increases with 
the number of photons involved. This indicates a direct 
emission process immediately induced by the pulse for 
high-order ATI. On the other hand, the more filled in- 
ner rings correspond to isotropic emission, which is most 
probably related to coUisional thcrmalization. A solid 
understanding of the PAD structures, however, remains 
a challenging task for future theoretical investigations. 

The above results indicate that electronic thcrmaliza- 
tion becomes increasingly important with increasing re- 
action time, i.e., if the laser pulse length exceeds the time 
required for coUisional relaxation. Indeed, the underlying 
energy distributions in Figs. pl| and 22 show a roughly 



exponential behavior. For Na,!^ perfect exponentials have 



been measured ( Bchlipper et al\ ^001|). The slopes s in 
exp (—siskin) of these distributions depend on the laser 
intensity, see the filled boxes in Fig. Clearly the slopes 
decrease with rising intensity, which can be interpreted 
as an increased heating of the system. As a matter of 
fact, the situation might be even more complex: TDLDA 
calculations (without electron-electron collisions) on the 
same Na^g also yield spectra with smooth exponentials 
while it is for certain that the signal stems from direct 
emission only. Trend and magnitude agree fairly well 
with the experimental data, cf. Fig. This indicates 
that an interpretation as thermal emission is not compul- 
sory from these data alone. A combined analysis includ- 
ing angular distributions and ideally also time resolved 
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FIG. 21 Measured photoelectron spectra from Ceo exposed 
to laser pulses (790 nm) for various pulse widths and inten- 
sities (as indicated). Dashed lines show e stimated thermal 
contributions. After (Campbell et al., 200C). 
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FIG. 22 Photoelectron angular distribution of Ceo for irra- threshold, cf. Sec 



diation with 800 nm laser pulses of intensity / = 10 W/cm 
and pulse width of 60 fs. The laser polarization is oriented 
along the horizontal axis. The rings correspond to above- 
threshold ionization, with increasing alignment for larger 
numbers v of absorbed photons. The inset gives the angle- 
integrated intensi ty, in qualitative agreeni ent with the result 
in Fig. ^ After ( |Skruszewicz et aZ.| , pOOsj ) 



VI. CLUSTER DYNAMICS IN STRONG FIELDS 

For the previously considered low and intermediate 
laser intensities, where nonlinearities like plasmon broad- 
ening, the onset of saturation effects, and ATI have been 
discussed, the electronic configuration and ionic struc- 
ture of the cluster remains in many respects at least sim- 
ilar to that of the initial state. This situation changes 
for laser-cluster interactions in the strong-field domain, 
see examples Fig. |l|c-d. In this regime the interaction 
leads to radical changes of the structure and proper- 
ties of the clusters, such as, e.g., light-induced metal- 
lization of rare-gas systems through strong inner ioniza- 
tion, the creation of multiple core- vacancies, or full clus- 
ter destruction via Coulomb explosion as a result of ex- 
treme charging. The cluster response is mostly (at least 
in the IR) field strength dominated, as may be deduced 
from the Keldysh pa rame ter or comparison with the BSI- 
[I.C| . However, such an assignment 
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FIG. 23 Slope of the PES from Na^g irradiated by a 200 fs 
laser pulse of 3.1 eV. Results are drawn versus intensity (as 
login (J) with J i n unit s of W/cm^). Results from TDLDA- 
MD ( jPohl et all |2004a|) are comp ared with the experimental 



results of (Schlippcr et ai, 2001) using comparable experi 



mental conditions. 



measurements would be required in order to unambigu- 
ously distinguish between direct and thermal processes. 
These effects are taken care of in the models particularly 
suite d for hig hly excited dynamics, i.e., VUU si mulation s 
(Sec. III.B.3| ), classical molecu lar dynamics (Sec. [II. B. 4), 
and rate-equations (Sec. III.C ). In the experiments, how- 
ever, a clear identification of thermalization effects re- 
mains difficult. Here imaging techniques as demonstrated 
above provide a promising tool for further investigations. 



based on atomic adiabaticity parameters may not be use- 
ful in all cases. The above phenomenological classifica- 
tion seems to better express what is widely understood 
under ,, strong-field cluster dynamics". 

Several illuminating early experiments in the nineties 
have sparked the rapid development of this topic. For 
an overview, we first concentrate on key phenomena and 
their possible relevance for technical applications. Subse- 
quently we review aspects of the underlying microscopic 
mechanisms according to the current status of knowl- 
edge. Along this line, the followi ng discussion is di- 
vided into two parts. Section VI. A highlights early sur- 



prises and experimental key results, ranging from mea- 
surements of energy absorption and emission of energetic 
particles to short wavelength radiation, and is restricted 
to a generic discussion of mechanisms and typical trends. 
Section VLB reviews routes towards a more detailed mi- 



croscopic understanding by closer relating theory and ex- 
periment and by pursuing more elaborate schemes like 
time- or angular resolved analysis. For a space-saving 
description of laser parameters, i.e., peak intensity /q, 
pulses duration t (FWHM ), and wavelength A, we use 
the compact notation (/q ; t ; A). 



A. Early surprises and basic trends 

1. Laser energy absorption 

A remarkable property of clusters in intense laser fields 
is very efficient energy absorption. At intensities of the 
order of 10^^ W/cm^ the average energy capture per atom 
can attain values of tens to hundreds of keV and by far 
exceeds that of atomic and molecular targets. Basically 
all of the violent processes discussed below have their 
starting point in this enhanced absorption, in conjunction 
with the absence of dissipation into surrounding material. 

The direct measurement of the absorption from the rel- 
ative loss of laser pulse energy in the interaction region 
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requires a high target density (10^'^ — 10"'^^ clusters/cm'^). 
This situation can be reahzed close to the nozzle (~1 mm) 
of supersonic gas expansion sources and corresponds to 
an effective particle spacing of >100nm. As a typical 
exarnple for the different behavior of clusters and gases, 
Fig. ^ displays the relative energy absorption Al of high 
intensity laser pulses in a dense jet of Ar^r compared to 
a Ne gas as function of backing pressure. Since Nc does 



worked out in Sec. VLB 




2. Highly charged atomic ions 

The strong optical absorption leads to high ioniza- 
tion and usually complete disintegration of the clus- 
ters. Atomic ions with high ionization stages q are fi- 
nally detected, see Fig. fo r an early ion spectrum 

After excitation 
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Gas Jet Backing Pressure (bar) 



on Xe clusters (Snyder et al. 

with (10^^ W/cm"^ ; 350 fs ; 624nm)-pulses a broad charge 
state distribution emerges, extending up to q™'^^ = 20. 
Such charge states are much higher than those from 
atomic Xe under similar conditions. For example, 
pulse intensities of 10^^ W/cm^ arc required to produce 
Xe^^+ (Dammasch et al, pOOl ) from atomic gases - in 
reasonable agreement with the BSI model. From 65 A Xe 
clusters, io n s with g™^^ = 40 were reported by ( Dit- 
pM7b| ) after exposure to laser pulses with 



mire et 



FIG. 24 Energy absorption of Arjv and Ne exposed to in- 
tense laser pulses (7 x 10^^ W/cm^ ;2ps;527nm) as a func- 
tion of backing pressure. The estimated Atm cluster diameter 
is 80 A at 40 bar, and 100 A at 55 bar. A major fraction of the 
laser pulse (up to 80%) is absorbed by the large Ar clusters 
while atomic Ne ga s rema ins nearly transparent. Adapted 
from (Ditmire et al, 1997E). 



(2 X 10^" W/cm^ ; 150fs ; TSOnm). The direct compari- 
son of Arjv (N ~ 100) to Ar gas after irradiation with 
(2 X 10 ^"^ W/cm^ ; 30ps ; 1 064nm)-pulses shows similar 



trends (Lezius et al., 1997), i.e., substantially higher max- 
imum charge states with clusters (^q"^^'^ — 10 with clus- 
ters over q™^'' — 3 with gas). The ion spectra from metal 
clusters (Ag^r, Auat, Ft at, Pb^r) were explored with fcm- 



not icondcnsc at room temperature, the overall low ab- machcr et al 



sorption in Nc provides a reference for a gas of similar 
atomic density. In contrast to that, Ar clusters become 
increasingly opaque with cluster size beyond the onset of 

ampin)? 



several studies 


KoUer et al. 


fladcliffc et al. 


2005: Schu- 



cluster formation (at about 5 bar in this particular ex- 
it rihould be noted, that laser pnlaa depletion 



1999| ), leading to values for g"^^^ up to 
30 for intensities below lO^^W/cm^. Further, highly 
charged ions have been reported for molecular clusters, 
i.e., a tomic iodine up to q = 15 from (CHs^at (Ford 



from light scatter ing was found to be insignificant ( Dit 



mire et al., 1997a). Further, the increase of gas flux with 
pressure can be ruled out as a major origin for the higher 
pulse depletion (Zweiback et al, 2002| ). Substantial ab- 



sorption of up to A™^^ = 0.8 for the highest pressure in 
Fig. ^ is typical for dense atomic or molecular cluster 
beams and has also been observed with (Il 2)jv, ( D2)jv, 



et al\ , |1999D with (2 x 10^^ W/cm^ ; 130 fs ; 795 nm), and 
up to 0''+ from to (Il20)Ar (Kumarappan et al, p0031: ) 
with (8 X 10^^ W/cm2 ; 100 fs ; 806nm)-pulses. 

As a general remark it should be noted that ion dis- 
tributions like the one in Fig. |2^ reflect an average over 
the focal intensity profile in the interaction zone where 



KrAf, Xcat (Ditmire et al, 1997a, 1999; Lin et al, |200]J 



Miura et al , 2001). In most cases, similar to Fig. |24 
laser attenuation scales roughly linear ly with st a gnati on 
pressure and then saturates, see also ( Jha et al , ^006 ). 

When analyzed as a function of pulse intensity, sub- 
stantial absorption sets in at relatively sharp thresholds. 
Beyond the onset intensities, e.g., Jth « 3 x 10^'^ W/cm^ 
for ArAT and Ith ~ 4 x 10^^ W/cm^ for Xcat at 527 nm, 
the pulse depletion increases rapidly and attains 
Ah ^ 0.5 A™'^^ at one order of magnitude higher intensity 



( Ditmire et al , 1997a). As the thresholds roughly follow 
the trend of the corresponding atomic BSI intensities, 
the behavior indicates an avalanche breakdown process 
triggered by atomic optical field ionization to establish 
efficient absorption. A clear signature of the dynami- 
cal nature of the energy capture, which turns out to be 
largely driven by resonant collective electron excitation, 
is the pulse length dependence ofAL- This will be further 
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FIG. 25 Mass spectrum of highly charged atomic 
Xe''^ ions resulting from by excitation of Xcat with 
(1 X 10^^ W/cm^ ; 350 fs ; 6 24 nm) showiiig char ge states up to 



— 20. Adapted from (Snyder et al, 1996) 
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regions of higher intensity contribute with a smaller ef- 
fective volume. Only recently, it has been demonstrated 
that contributions from the different intensities to such 



of Sec. VII. B. For now we come back to excitations with 



spect ra can be deconvoluted ( Doppncr et al , 2007b 



2009) by using intensity-selective scanning. 

Dedicated studies on heterogeneous, doped, and 
embedded clusters have been performed to in vesti- 
gat e the ef fe ct of the cluster composition. In ( Pur- 
ncll et al.. 1994), irradiation of HI clusters 



al. , 1994), irradiation of HI clusters with 
(1 x' lO^^' W/cl?T350 fs ; 624 nm) yields I«+ with up to 
q = 17. In the same work, Ar^v and Ar atoms attached 
to (HI) AT yield Ar''+ up to g = 8 for Arjv(HI)M, whereas 
no notable contribution from multicharged Ar''"'" ions 
is found for bare ArAr. This supports, that the low- 
IP atoms or molecules act as chromophorcs and initiate 
nanoplasma formation. Subsequently also constituents 
with more strongly bound electrons can be ionized, e.g., 
via electron impact ionization. Other matrix effects oc- 
cur in helium nanodroplets: Experiments on embedded 
clusters have shown evidence for electron transfer pro- 
cesses, where highly charged ions capture electrons from 
the surrounding helium (Doppncr et al., 2007a). When 



considerably ionized, the helium shell can produce strong 
absorption e nhancement due to reson ant heating of the 
nanomatrix ( Mikaberidzc et al, 2008). 

Whereas most of the results have been obtained with 
optical lasers, first experiments are at hand making 



use of a VUV Free Electr on Laser ( [Laarmann et al. 
2004|; IWabnitz et al\, |2002[). Power densities of up to 



3 X 10 ■^W/cm-^ at 98 nm (=12.65 eV) were used in these 
experiments on rare gas clusters. Note that IBS heating 
is less effective at shorter wavelengths because of the low 
ponderomotive potential so that multip hoton ionization 
conditions are expected, cf. Sec. II. C| . Moreover, reso- 
nant collective heating can be ruled out because of the 
high laser frequency. Still, ionization of clusters is quite 
effective, leading to ions with charge states up to Xe^"'" 
and Ar^"'". These findings underline that strong cluster 
excitation is still possible in the domain of large Keldysh 
parameters. The origin of the high energy absorption 
required for the observed charging has been investigated 
by several groups. Various concepts were proposed, rang- 
ing from models based on enhanced IBS heating due to 
strong electron-ion scattering (Bantra and Greene, p003 ), 
over efficient IBS heating resulting from a high-density 
nanoplasma produced by local fiel d enhancement of in- 
ner io nization by neighboring ions ( Sicdschlag and Rost , 
2004), to many body heat ing effects ( Bauer , |2004b| ; Jun- 
grcrthmayer et al, 2005). It should be noted that the 



calculations performed by ( piedschlag and Rost| , |200J) 
for Xcgo show good agreement with the experimental ion 
spectra when including the experimental focus averag- 
ing. Within the approach of (Santra and Greene, 2003| ) 
there remain clear deviations from the experimental ion 
spectrum when taking the focus effect into account. For 
more details see a recent review with emphasize on the 

p006| ). Further perspec- 



VUV-domain (Saalmann et al 



optical lasers. 



3. Ion energy distributions 

Another early surprise was the large kinetic energy 
of atomic species emitted from clusters in intense laser 
pulses. Due to the strong heating of cluster electrons 
and high cluster charging on the femtosecond time scale, 
huge amounts of thermal and Coulomb energy are avail- 
able to be released within in the explosion of the system. 
In experiments atomic ions from Xcjy with ki netic ener- 



gies beyond 1 MeV were observed ( Ditmirc et al , 1997c). 
This has opened the route to tabl e -top e xperiments on 
cluster-based fusion (Ditmirc et al, 1999). Interestingly, 
for rare-gas clusters a rather sharp onset of high energy 
ion emission is observed. For XeA? a thresh old intensity 
some what above lO^'^W/ cm^ was reported ( [lisch et al 



2003) with 230 fs pulses at 790 nm being roughly compa- 
rable with the BSI threshold intensity, see Eq. (^. 

The charge-state averaged ion energy distributions 
turn out to be very broad, see Fig. Wd for Pb^r and 
Fig. ^ for (N2)Ar, and show a smooth decrease with 
increasing energy, often followed by a cutoff which is 
frequently termed "knee" -feature. Typically the maxi- 
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FIG. 26 Measured ion energy spectra from {N2)n (curves ) for 
excitation with (1 x lO^'' W/cm^ ; 100 fs ; 800 nm)-pulses ( |Kr^ 
ishnamurthy et al, 200^ ). The highest stagnation pressure 
corresponds to (A*') — 2300. With increasing backing pres- 
sure, i.e., for larger clusters, the spectra are shifted to higher 
energies. Symbols represent fits using a Coulomb explosion 
model that incorporates averaging due to t he laser beam pro - 



file and the cluster size distribution. After ([slam et al, 2006) 



fives of VUV and XUV excitations of clusters are subject 



mum energy, which may be quantified by the energy of 
the knee-feature, increases wit h cluster size, se e Fig. 2^ . 
However, in the case of Xcat (Mendham et al, 2001) it 
was found that the maximum ion energy grows with clus- 
ter size until it levels out. For sufficiently short laser 
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pulses the nanoplasma model predicts that the ion energy 
decreases beyond a certain size, since the slower expan- 
sion ofjajge_clustcrsJm£cd£Sj^ 



of a considerable fraction of highly charged heavy- 
element ions in the cluster pro duces a strong l y repelling 
background for t h e light i ons (iGriUon et ali ^002| ; ^ 



The saturation can be explained by the relatively broad marappan et al., 2003b; Madison et al. , 2004). For 



experimental cluster size distribution, i.e., the molecular 
beam still contains optimally sized particles producing 
the maximum ion energy. 

Several additional effects contribute to the shape of the 
ion energy spectra, i.e., the spatial laser intensity profile 
and the degree of cluster ionization. Taking all these ef- 
into account, experimental data on Xcjy (|Ditmire 



feet! 
et a. 
et a. 



ishnamurthy et al. 



1997c| ; ppringate et al. , 2000f|), Ar jy ( pK^umarappan 
20011), (H2)jv (ISakabe et alj |2004|), and (N2)2 (|K~ 
lool " can b e reasonabl y well fitted 
considering Coulomb explosion ( Islam et al. , 2006 ). 

Neglecting thermal electron excitation and assuming a 
uniformly charged monoatomic cluster, the final kinetic 
energy e/ of an atomic ion is de t ermined by its initial 



potential energy ( Last et al. , 1997 ; Nishihara et al. , 2001 
Zwciback etd\ ^OOOj) 



An 



2 2 

pir q 



14.4eVA, 



(22) 



where r is the initial radial ion position, q the charge 
state, and pi the number density of ions in the cluster. 
Hence, ions at the cluster surface acquire the highest 
recoil energies e™^^ and this maximum energy also in- 
creases with cluster size. The latter trend was experimen- 
tally confirmed by several groups, e.g., on Xejv and Ar^v 



by (iDitmire et al\, |l997b|; [Lezius et aLj |1998| ; |Li et al 
2003bD . In ( IKrishnamurthy et q/.| , |2004| ), a monotonous 



rise of ef'^^ from about 1 keV to about 



.5 keV was found 
with (N2)Ar when increasing the size from A'^ = 50 to 
N = 2300, cf. Fig. |26[ With Pb^v, ef'''' rises from 70 to 
180 keV whe n increasing t he clu ster size from A^ w 100 
to N 500 (iTcuber et alj ^OOl]) , see Fig. |l|d. 



At constant cluster size, a recoil energy enhancement 
is observe d when ad d ing spurious amount s of ap propriate 
dopants ( |Jha et al\ , ^006| ; [PurneU et al] , |1994| ). For in- 
stance, Ar clusters (N ~ 2000) containing abou t 60 H2O 
molecules were considered in (|Jha et al. , 2006). Under 
exposure to pulses with (1 x 10 W/cm"^ ; 100 fs ; 800 nm) 
the ion yield at 100 keV scales up by a factor of three 
and the maximum ion energy is larger when compared 
to the dopant-free case. This effect was traced back to 
a longer phase of strong cluster heating and ionization, 
since dopants with lower ionization thresholds can initi- 
ate nanoplasma formation earlier in the laser pulse. 

Clusters containing a mixture of low and high atomic 
number elements can be used to enhance the kinetic en- 
ergy of the low atomic number ions. This is of par- 
ticular interest for the acceleration of H+ or D+ for 
fusion reactions. For pure (D2)Ar (A^ < 10^ ) ion ki- 
netic energ ies up to 30keV were detected (Zwciback 
et ^002| ) with (1 x 10^^ W/cm^ ; 35 fs ; 820 nm)-pulses 
and energies up to 8.1keV we r e foimd for bare 
(H2)Ar {N ~ 10^) dSakabe et |2006i ^004] ) with 
(6 X 10^6 W/cm2 ; 130 fs ; 850 nm)-pulscs. The presence 



(D20)Ar an enhancement in e/(Z3+) of 5.6 over the re- 
sult from (D2)7v of the same radius is predicted in the 
limit of complete and instantaneous ionization due to 



higher ionization stages of oxygen (Last and Jortner 
2001). The increase in the kinetic energy of D_J__was 
verifie d in a study on (D2)Ar and {CTli)^ (Madison 



et al., 20041 ). From deuterated methane excited with 
(1 X 10^ ' W/cm^ ; 35fs ; 820nm)-pulscs, deuterium ener- 
gies of up to 120 keV were found (Grillon et al., 2002). 
Interestingly, only doubly charged carbon was detected, 
indicating substantial electron recapture. This is fur- 
ther supported by the fact that the maximum C"*" energy 
(180 keV) substantially exceeds the value for D+. 

Additional insight into the explosion dynamics and the 
initial ion position within the cluster can be gained by 
simultaneously measuring charge states and energies. In 
principle, each emitted ion state has its own character- 
istic spectrum. To access the charge-resolved spectrum, 
techniques such as MD-TOF or Thomson spectroscopy 



can be applied, see Sec. fV.C. As an example obtained 



with another method, i.e., retarding field analysis. Fig. |27| 
presents charge-resolved spectra measured at various re- 
coil energies e/ from irradiation of Xe2.5oo with pulses 
of (2 X 10^6 W/cm2 ; 150 fs ; 800 nm). With increasing e/, 
the charge state distributions shift and broaden. At 
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explosion 



Charge-resolved ion 
of Xe2500 exposed 



spectra from Coulomb 
to laser pulses with 
(2 X 10^'^W/cm^;150fs;780nm) for different recoil en- 
ergies ti (as indicated). With increasing ion kinetic energy 
the charge sp ectrum shifts to hi gher ionization stages. 
Adapted from (Ditmire et al., 1997b). 
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100 keV, ions with q = 2A are the most numerous and 
charge states up to g ~ 40 are observed. In stu dies on 
Ar^v [N = 1.8 X 105) and Xe^ = 2 x 10*^) ( |Lezius 
et q||, |1998D a scahng Uke e/(g) = 180 eVg^ (Ar) and 
ei{q) = 160 eVq^ (Xe) was found for q <6, as expected 



from electrostatic consideration, see Eq. < \22\) . The higher 
charge states (q > 10) show a more hnca r dependence on 
q. A similar behavior was reported by (Lebeault et al. 



2002). This linear dependence has frequently been in- 
terpreted as a distinct fingerprint from hydrodynamic 
cluster expansion (driven by thermal electron energy), as 
predicted by the nanoplasma model, see Sec. III.C. Such 
an assignment of parts of the spectrum to a Coulomb 
explosion or a hydrodynamic expansion has neverthe- 
less to be made with care. From a theoretical point of 



view (Ditmirc et al., 1996), the conclusion that hydro- 



dynamic forces dominate the expansion is based on the 
assumption of a constant electron temperature in the ex- 
panding nanoplasma. This might be too crude to de- 
scribe the dynamics correctly, as the nanoplasma experi- 
ences efficient expansion cooling during cluster explosion, 
see Sec. VI.B.l. A more elaborate analysis of the expan- 
sion process predicts a scaling law for the maximum recoil 
energy as a function of atomic density, c luster radius and 
initial temperature (Pcano et al., 2007, 2006| ). A one-to- 
one assignment of a quadratic dependence of (q) with a 
Coulomb explosion, see Eq. ( p^ ) , is correct only if the ion- 
ization process is quasi-instantaneous, i.e., much shorter 



than t he tim e scale of the ionic motion, see also (Tcuber 
et |200l| ). For a dynamical cluster charging during 
the expansion, even pure Coulomb explosion can lead to 
a linear ei{q) scaling. Note that this interplay between 
charging time and ionic motion has also been discussed in 
the context of structure analysis of biomoleculcs through 
scattering of XFEL radiation (Neutze et al., 2000). 
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FIG. 28 Top: K-shell line emission from Axn {{N) ~ 4 x 10*') 
exposed to ( 1.6 x 10^^ W/cnf ; 500f s ; 800 nm) laser 

Bottom: L- 
excited with 



pulses. From (Dorchies et al., 2005) 



shell line emission from Xejv {{N) — 12) 



lO^iW/cm^ ; < Ips ; 800 nm/248 nm). From ( Schroedei 



et al 



1998), with kind permission from lOP Publishing. 



4. Soft x-ray and EUV emission 

A decay channel of relevance for diagnostics and possi- 
ble applications of laser-cluster-interactions is the emis- 
sion of energetic photons, particularly in the soft x- 



ray idomalii. — The x-ray specLra coiiLalii Hue ciiilHsloii 
that' iLiflecLs lecunibhiaLiuu uf ulecLiuiis iii weakly buLuid 
atomic levels with core-level vacancies, as first r e porte d 
by Rhodes and co-workers ( McPhcrson et al. , 1994). 
Note that the charge distribution during the pulse, as 
partly reficctcd by the x-ray spectra, may be different 
from the final ion spectra because of free-bound electron- 
ion recombinations and charge transfer processes. 

The two examples in Fig. ^ display such line emis- 
sion spectra that can be attributed to specific ion charge 
states. The to p panel gives a result on K-shell emis- 
sion from Ar^r (Dorchies et al, 2005) for laser intensi- 
ties of lO^^W/cm"^. The example in the lower panel has 
been obtained at 10^^ W/cm^ on small Xejv and exhibits 
x-rays down to 2.4 A(5.2keV) a nd contributions from 
highly charged Xe«+ up to q 40 ([Schroeder et a/.|, |2001 



1998| ). Note that the excitation at 248 nm compared to 
800 nm results in lines from much higher charge states. 
So far, the most energetic x-ray emission from nano-sized 
targets was observed on Krjv where strong i^a,fj radiation 
(12.66keV/ 1.02A and 14.1 keV/0.88 A) was found ( psac 



et al. 



2004), sec also Fig. |i|c. A comparative study on 
Kr AT and Xcat with 5 x 10^^ W/cm^ concentrated on 



Ar jv ; rvr^r ana a.c n 
M-shcll line emission, e.g., Kr^ + (3d 



-4d), in th e wave - 



length regime of about 10 nm ( McPherson et al , 1993). 
Actually this was the first evidence that ionization dy- 
namics in clusters under strong field conditions signifi- 
cantly differs from the behavior of atomic targets. The 
enhanced x-ray emission disappears when the clusters are 
destroyed by a weaker pre-pulse ( ^kobelev et al , 2002). 
The major contribution of the x-ray photons is emitted 



on the ns time scale (Ditmirc et al, 199£; Kondo et al. 



2002; Larsson and Sjogren , 199E ), which is comparable 
to recombination lifetimes. 

The yield of x-ray emission with laser intensity scales 

as 

/3/2 

once the saturation regime is attained. This be- 
havior reflects the effective volume in the focus (Dobosz 
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Laser peak intensity [W/cm 



FIG. 29 Absolute x-ray yield versus the laser intensity Xe]v 
{{N) ~ 10^) for two different pulse durations. Solid lines re- 
flect the effective focal volumes with intensity thresholds of 
3.5 X f OiiW/cmifor 60 fs and 2.5 x f O^"* W/cm^ for 320 fs. 
From (Lamour et al, 2007), with permission from lOP Pub- 
lishing. 



et at, 1997 ; Rozct et ai , t-'nn|| i and is similar to atomic 



targlets ( jAugustc et al\ , |1992D : On KrAr (N~ iO'') a 
threshold for ihe onsci of high energy phoron emission 



was observed slight ly below 10 ° W/cm- 



130 fs pulse width (Dobosz et al, 1997) 



at 790 nm and 
The onset in- 
tensities, however, show a strong dependence on pulse 
duration and cluster material. Pulse length dependent 
and X cat were ob served in (Lamoui 



thres holds f or Ar jy 

et a|], |2005| , ^007| ; [Prigent e/: al] , [20081) . An example from 



Xcat is displayed in Fig. |29|, where threshold intensities 
of 3.5 X lO^^W/cm^ are found with 60 fs pulses, whereas 
320 fs pulses result in a substantially lower threshold of 
only 2.5 x 10^''W/cm^. T he data hint at a transiently 
resonant heating, see Sec. VI.B.l . Note that the low 



threshold for the long pulse case is only slightly higher 
than the BSI threshold for atomic Xe (9 x 10" W/cm^). 

A high number of emitted photons is crucial when 
aiming at technical applications like EUV lithography 
(EUVL), see e.g. ([Banine and Moors], |2004|). Further- 



more, for EUVL applications the emission must occur in 
a narrow spectral range to avoid aberrations. The spec- 
tral range around 13 nm is of importance since multilayer 
mirrors of M o:Be and Mo:Si re ach high reflectivities of 
nearly 70% ( ^tuik et 



19991 ). Conversion efficiencies 



of about 1%/ 27r sr have been achieved at only 2 % band- 
width, underlining that line emission fr om laser-driven 



Xe clusters, droplets, or jets, see e.g. (Hansson et al. 



2004), might by suita ble for ne x t generation micropro- 
cesso r manufacturing (Attwood, 2007 ; Wu and Kumai , 
2007). Within the soft x-ray regime, a recent study has 
demonstrated the high potential of clusters in intense 
laser fields as debris-free radiation sources for nanostruc- 



ture imaging ( Fukuda et al , 2008 ) 




High liarmonic generation 

High Harmonic Generation (HHG) has been q uite ex- 
tensively studied in atomic and m olecular gas jets ( Chang 



1997 



K rausc et al 

pool 



1992| ; ppiclmann et al\ |1997| ; 



Vclotta et al\ p001| ). The resulting spectra contain odd 
harmonics because of inversion symmetry and show an 
initial strong intensity decrease, a plateau-region, and 



a rap id cutoff near £'ip-|-3.17Up (Brabec and Krausz 
2000| ). This cut-off is governed by the maximum return 



energy of electrons ( porkuni , 1993) and reflects the im 
portance of coherent st imulated recombina tion in atomic 
and molecular targets ( [Pukhov et al , 200S ) . The physics 
of HHG in clusters and particles has been investigated 
only in a few experiments, although clusters may act as 
a unique nonlinear optical medium. 

In an early study on HHG in clusters, a substantial 
harmonic signal of high order, actually 23th harmoni c 
(Htf'"'), was reported for ArAr ( [Donnelly et al[ [1996 ). 
For laser intensities of up to 1.5 x 10^* W/cm-^, the HH'^'^ 
signal scales with l'^''^^ and then changes its slope to 
I"**^. In atomic gases for comparison, a po wer depen- 
dence of I^^ is found in the cut-off region (Wahlstrom 
et ai , 1993). A clear enhancement effect of the HHG 



yield in clusters was also reported on Xcat (Tisch et al 



1997), where, for HH^ (78 nm), an increase of almost one 



order of magnitude was found with respect to the atomic 
gas. A comparison of HHG-spectra from atomic argon 
and Ar^v taken under identical laser conditions is shown 
in Fig. Close to the nozzle, (i.e. a; = 0, see inset), the 
cluster formation is not completed and HHG mainly re- 
sults from atoms. When clusters are present in the beam 
(at X = 1.5mm), the cut-off wavelength decreases consid- 
erably from 17 nm to 14 nm, thus increasing the highest 



xk 1.5 mm 




30 27 24 21 18 15 

Wavelength (nm) 



FIG. 30 High harmonic signals from Ar vs. Arjv: spec- 
tra resulting mainly from atomic argon (dashed curve) 
and Arjv ((.^V) ~ 10^, full curve) exposed to pulses with 
(2 X 10" W/cm^ ; 25 fs ; 800 nm). The probed target condi- 
tions are selected by shifting the laser focus with resp ect to 
thenozzle (as indicated). Reprinted with permission (Vozzi 
et al, 2005). Copyright 2005, American Institute of Physics. 
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HH order beyond the ponderomo t ive li mit. A similar re- 
sult was obtained in (Pai et ai, 2006). Common with 



all results is a significant HHG enhancement in clusters 
when compared to an atomic gas. 

In fact, the special linear wave propagation properties 
of dense cluster media (cluster-cluster separations of the 
order of the diameter) could be advantageous for efficient 
HHG. Since a gas of inner-ionized clusters can build up 
strong depolarization fields, electromagnetic waves can 
propagate below the plasma cutoff in a particular optical 
mode-the so-called cluster mode (Tajima et al, 1999). 
In a homogenous plasma, waves with frequencies below 
the plasma frequency become evanescent. Thus, a much 
higher electron density can be established in the cluster 
media without reflecting the fundamental wave. More- 
over, in contrast to atomic gas plasmas, the refractive 
index of cluster media can be larger than one. A cer- 
tain mixture of atoms and clusters can be used to tune 
the refractive index in order to fulfill a major require- 
ment for efficient HHG, i.e., phase-matching. In princi- 
ple phase-matching can then be attained for any desired 



harmonic (Tisch, 2000), underlining the promising pos- 
sibilities of clusters for tailored optical media. 

Nonetheless, details on the microscopic mechanisms of 
HHG in clusters, e.g., concerning the interplay of stim- 
ulated recombination and bremsstrahlung, are yet to be 
explored. For a theoreti cal study on the contrib ution 



from bremsstrahlung see ( [Popruzhenko et al\ , p008| ) 



B. Analyzing the microscopic cluster response 

The above examples highlight the violent and multi- 
faceted nature of laser-cluster interactions in the strong- 
field regime. The aim of the following sections is to review 
selected aspects of the ultrafast microscopic dynamics for 
excitation with strong optical lasers in more detail. This 
concerns the dynamics of cluster heating, expansion, and 
ionization. The influence of the pulse structure on the 
laser energy absorption as well as on the emission of elec- 
trons, ions, and x-rays has be studied in great detail in 
experiments with stretched pulses or by dual-pulse exci- 
tation. In most cases the observed signatures support a 
large impact of resonant plasmon excitations as will be 
discussed in connection with theoretical concepts and re- 
sults from numerical simulations. Further, we comment 
on the persisting difficulties in explaining the origin of 
highly charged ions, which requires to understand the 
strong-field ionization on the atomic scale in the pres- 
ence of a highly excited many-body environment. Fur- 
ther, angular resolved emission of electrons and ions will 
be addressed, which reveals unique acceleration effects 
in laser-excited clusters. Such studies can help to iden- 
tify sensitive parameters for a control of specific decay 
channels and are important tests for theoretical models. 
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FIG. 31 Laser power absorption (symbols) by Xe clusters of 
different size (as indicated) for dual-pulse excitation (a) and 
irradiation with stretclied pulses (b) at A=810nm. The dual- 
pulse results were obtained with 50 fs pulses of peak intensities 
/i = 1.6xlO"W/cm2 and72 = 1.6 x lO^^W/cm^ The single- 
pulse experiment was performed with constant pulse energy 
of 6.5 mj, resulting in a peak intensity of 2.3 x lO^^W/cm^ 
at 50 fs pulse duration. The curves show results calcu- 
lated with the hydrodynarn ic model from Ditmire. Adapted 
from (Zweiback et ai, 1999). 



1. The l<ey role of collective excitations 

a. Evidence for resonant absorption: An important 

concept for explaining the high energy absorption of clus- 
ters in intense optical laser pulses involves resonant col- 
lective electron excitations. Since an inner ionized clus- 
ter is usually overcritical, a certain density lowering from 
cluster expansion is required to achieve frequency match- 
ing of th e collective mode with the IR laser pulse, see 
Sec. II. C. Striking experimen tal evidence for such tran- 



19991) , see 



sient resonance was found in ( Zweiback et al 
Fig. where the laser energy absorption in a Xe clus- 
ter beam is measured as a function of pulse profile. The 
upper panel shows the result of a dual-pulse (pump-and- 
probe) experiment. In a simplified picture, the leading 
pulse excites the cluster moderately and initiates its ex- 
pansion. At a certain time the system reaches resonant 
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conditions leading to a strong peak in the absorption as 
is clearly seen from the figure. The optimal pulse delay 
increases with size, showing that larger clusters require 
more time to reach frequency matching. An alternative 
way to explore the expansion time is to use one pulse but 
varying its length, as demonstrated in Fig. pl|b. Simi- 
larly, it shows an optimal pulse duration to induce maxi- 
mum absorption which again increases with system size. 
A word of caution is advised for the interpretation of 
such stretched pulse measurements, as the intensity de- 
creases with increasing pulse duration. For example, con- 
sidering a certain threshold intensity for the nanoplasma 
buildup (e.g. the BSI intensity), the effective interac- 
tion volume is strongly dependent on the pulse duration 
as well. Thus, a long pulse with low peak intensity will 
probe a smaller number of targets. This problem docs 
not occur with a dual-pulse setup as used in the case of 
Fig. pT|a. Nevertheless, both excitation schemes show a 
pronounced resonance feature in qualitative accordance 
with calculations based o n the h ydrodynamic model from 
Ditmire discussed in Sec. III.C, see curves in Fig. |3^. 

Modeling the cluster response by a single collective 
mode is, however, strongly oversimplified. A more com- 
prehensive picture can be dra wn from microscopic simu- 
lations, such as MD, see Sec. [II. B. 4. As an illustrative 
example. Fig. ^ displays results from a MD simulation 
of Xe5083 exposed to a 250 fs pulse of lO ^^W/cm^, for 



„ 1.0 



methodic details see ( [Fennel et al.[ |2007bD 



The upper panel shows the effective absorption cross- 
section (solid line), as derived from the total energy cap- 
ture (dashed line), together with the laser intensity enve- 
lope (centered at time zero) . The broad absorption peak 
between —80 and 50 fs corresponds to the expansion- 
induced collective resonance. The dominant energy cap- 
ture proceeds near the crossing of the mean cluster radius 
R with the critical radius for resonant collective coupling 
i?crit at t w Ofs, see panel (b). The value of i?crit is 
estimated from the Mie formula (Eq. [^) by 

X 1/3 

i?crit - ( -{q)NX'' , (23) 



with (q) the average inner ionization of ions. For exam- 
ple, at A = 800 nm, a xenon cluster (rg = 2.5 A) must ex- 
pand by a factor of 1.4 for (q) = 1 and of 3.0 for {q) = 10 
to become resonant. The evolution of cluster inner ion- 
ization, see panel (d), leads to a time dependent Rcrit- 
In particular, rapid inner ionization triggered by optical 
field effects and further enhanced by electron impact in- 
duces a sudden rise of Rcrit in the leading edge of the 
pulse, see panel (b). That is accompanied by a short 
period of resonant coupling when -Rcrit crosses the ac- 
tual system radius before the cluster becomes overcrit- 
ical. This early ionization-driven resonance is reflected 
in a small feature in the cross-section at t Ki —240 fs and 
was a lso ob served in MD calculations by ( |Saalmann and 
Ros fl, |2003D. Such early resonance would not occur in 
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met il clusters probed with IR pulses, as they are over- 
critical already in the initial state. However, the major 



FIG. 32 Simulated dynamics of Xesoss exposed to a laser 
pulse with dO^^ W/cm^ :250 fs ; 800 nm) using the MD code 
from ( Fennel et al. , 2007b ). (a) Pulse intensity profile, 
cycle-averaged energy capture Sabs, and associated absorp- 
tion cross-section a^n — E^hs/I- (b) Mean cluster radius 
R — -Rrms, where i?rms is the root-mean-square radius, 

radial position of outermost ion Roi , and critical radius i?crit 
for resonant coupling, Eq. (|2^). The vertical dotted lines mark 
the matching between R and 7?crit, which coincides with the 
maxima in the absorption cross-section. The difference be- 
tween R and Roi indicates inhomogeneous cluster expansion, 
(c) Average and maximum electron kinetic within the cluster 
radius R. The dotted line shows a fit for adiabatic expansion 
cooling (see text), (d) Number of inner- ionized and contin- 
uum electrons (cycle-averaged); the difference of these values 
reflects the number of quasifree electrons, cf. Fig. ^ 



contribution to the energy capture i^abs proceeds within 
the expansion-driven resonance around tw fs. In total, 
the ultimately absorbe d energy exceeds 50 keV per atom, 
similar to findings of (BaalmanE, 2006). Note that the 

maximum cross-section of 0.75 per atom is compara- 
ble to values typical for collective resonances in the linear 
regime, cf. Fig. |l|b. The cross-section in Fig. ^2| a drops 
quickly after the resonance, reflecting the suppressed cou- 
pling efficiency at undercritical density. Collective effects 
become unimportant and IBS-heating of residual elec- 
trons is weak due to rare electron-ion collisions. 

Several effects contribute to the observed width of the 
resonance. To some extent, the broadening ca n be linked 
to inhomoge neous cluster expansion, sec e.g. (Milchbcrg 



et al. , 2001), with time-delayed resonant absorption in 
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radial shells of critical density. A comparison of the ra- 
dial position of the outermost ion i?oi with R in Fig. ^2| b 
show s that outer ions indeed exp and more quickly, see 
also ([ishikawa and Blenski , 2000 ). This is due to a less 
effecti ve scr eening; of ions near the cluster surface ( Peano 
et a Besides the influence of the ionic density 

profile, the driving of electrons beyond the cluster sur- 
face introduces an additional broadening due to nonlin- 
ear d amping ( jJungreuthmayer et al , 2004; Megi et al. 



2003| ) . For sufficiently high laser intensity, also the occur- 



ren ce a nonlinear resonance has been discussed (Kundu 
and |Bauei| , |2006t |Mulser et al\ , ^00^ . 

Coming back to the example in Fig. |3^, the average 
kinetic energy of cluster electrons of up to 1 keV (dashed 
line in panel (c)) shows that there is strong thermal exci- 
tation of the nanoplasma near the resonance. The max- 
imum electron kinetic energy within the cluster of up 
to 6keV (solid line) provides a reasonable measure for 
the depth of the transient space-charge potential pro- 
duced hran_oute£j£niMik)naj^ thermal excitation, see 
also ( ^aalmann and Rost , 2005). By electron impact 
excitation, such energetic electrons can directly create 
deep inner-shell vacancies required for hard x-ray emis- 
sion. After the resonant heating the expansion of the 
ionic background leads to an efficient electron cooling. 
This process can be well described by adiabatic cool- 
ing of an ideal gas in an expanding spherical vessel: 
(ekin(i)) = aR-^{t) + b, see dotted line in Fig. He. The 
offset parameter b accounts for the kinetic energy of elec- 
trons that become localized in ionic cells during expan- 



sion. As displayed in the lowermost panel, a substantial 



fraction of the inner-ionized electrons cannot be accel- 
erated to continuum energies and remains bound in the 
cluster potential. The further evolution of these electrons 
is of central importance for predicting the final ion charge 
spe ctrum and will be reconsidered in Sec. Vf .B.2|. 



T |ie above MP analysis illustrates that the dynamics 
is clearly dominated by collective energy absorption near 
the critical density. The cluster expansion to Hcrit there- 
fore sets a crucial time scale for strong-field laser-cluster 
interactions in the IR regime. Only at very high inten- 
sities, where the laser field exceeds the restoring force 
from the ion background potential, resonance effects can 
be disregarded (Heidenreich et al., 2007; Krainov and 
Smit-nov|, |2002|). 



b. Signatures in emission spectra: Having identified 
the dominant role of collective excitations in the absorp- 
tion, we now concentrate on corresponding signatures in 
the emission of x-rays, highly charged ions, and electrons. 
Time-resolved measurements of the x-ray emission 



stretched pulses ( 


Chen et al. 


, 2002< 


5; 


[ssac et al. 
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2004 


Lamour et al., ^005; 


Parra et al. 




2000 


; Prigcnt et al., 


2008; Zweiback et al. 


I1999p. 


An example for Xe clusters 



is given in Fig. |33| displaying the x-ray yield from 3d 2p 
transitions of Xe'^>^'*+ (Lamour et al., 2005| ). The signal 



can be interpreted as a measure of energetic cluster elec- 
trons, as the production of Xe^''"'" plus the 2p- vacancy by 
electron impact requires considerable energies of 7.3 keV 
and ^4.5 keV, respectively. The energy for the creation of 
the vacancy should be transferred within one single col- 
lision event. Fig.^ shows a steep increase in the x-ray 
yield for pulse durations up to 250 fs, indicating a grow- 
ing number of multi-keV cluster electrons. Note, that 
this is compatible with the generation of keV electrons at 
the instant of resonant heating in Fig. ^2| c. The optimal 
duration thus indicates efficient collective heating, see 



also (Parra et al., 2000; Zweiback et al., 1999) and similar 



experiments on EUV emission ( phen et Q^X|2002a| ). An 



alternative electron heating mechanism, namely multi- 
ple large-angle electron-ion backsc attering in phase w ith 



the laser field, was proposed in ( peiss et al. , 2006) 
order to explain the x-ray production from pulses that 
are too short for reaching resonant conditions. Whereas 
basic aspects of short wavelength emission from clusters 
can be rationalized, there are still several pending ques- 
tions. For example, the role of ionization and excitation 
over intermediate states or the impact of multielectron 
collisions on the production of core vacancies have not 
been resolved so far. Therefore, the physics behind x-ray 
emission from clusters remains a fascinating subject for 
further studies. 

Measurements of ion kinetic energy spectra as func- 
tion of pulse duration substantiate the strong impact of 
the pulse structure. Fig.^ shows a result of a constant 
peak intensity experiment on Xcat performed by ( Fukuda 
2003| ) . The mean ion energy grows with pulse dura- 



et al. 



tion reaching a maximum at about 500 fs, in accordance 
with the picture of a delayed resonance. A similar behav- 
ior has been reported from a constant fiuence measure- 
ment (Kumarappan et al, 2002). Also a clear effe ct of the 
temporal ph ase of the pulse has been observed (Fukuda 
2003). Significant differences in the ion energies 



et al 



were found for positively and negatively chirped pulses. 
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FIG. 33 X-ray yield at 4.4 keV from Xejv with (W) 4 x 10^ 
irradiated with stretched laser pulses of fixed energy (35mJ) 
at 800 nm wavelength. The shortest pulse (r=50fs) corre- 
sponds to a peak intensity of 3 x lO^^W/cm^ . The dotted 



line is guide to the eye. After (Lamour et al, 2005), Copy 



right 2005, with permission from Elsevier. 
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FIG. 34 Mean recoil energy of atomic ions emitted from Xejv 
clusters ((A'') ~ 5.5 x 10*) for excitation with stretched pulses 
at 800 nm (spectral width ^ 60 nm) and constant peak inten- 
sity of 2 X lO^^W/cm^. The results have been obtained with 
positively and negatively chirped pulses (as indicated). Note 
that a negative chirp corresponds to a decreasing laser fre- 
quency within the pulse. For 50 fs pulses the chirp rate is 
about 0.13nm/fs. Adapted from (Fukuda et al, 2003). 
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For negative chirp, i.e., a decreasing laser frequency with 
time, a 60% enhancement of the mean ion energy was 
observed, see Fig.|3^. This effect can be explained quali- 
tatively by the joint gradual frequency red-shift of both 
laser pulse and resonance, which, in turn, extends the 
time span for resonant collective absorption. Applying 
the nanoplasma model to these particular experimental 
parameters, 500 fs pulses with negative chirp lead to 1.2 
times higher total energy capture than the corresponding 
result with positive chirp. 

Evidence for an enhanced cluster ionization for cer- 
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FIG. 35 Charge state distribution of atomic ions emitted 
from lead clusters after exposure to laser pulses of variable 
duration and constant energy (19 mj). The laser peak inten- 
sity is 2.6 X lO^^W/cm^ for the shortest pulse (150 fs) and 
lO^^W/cm^ for the optimal duration (800 fs), where the lat- 
ter yie lds the highest c harge states of up to 28. Reprinted 

with kind permission from 
ledia. 



from (Doppner et al 



e.g. ( 


Doppner et al., 


2007b, 


2000; 


J^ukuda et al. 




20031 


Kollcr et al. 


, 199C; 


Lcbcault et al 


, 2002; 


Schumacher 


et a\.[ 1999). As a typical result. Fig. |35| displays spectra 



of high-q ions as a function of pul se width for an experi- 
ment with constant laser fluence (Doppner et al, ^000 ) 



Thclshortcst and most mtcnse pulses (lt)U tsj yield atomic 
ions up to q — 'zu. With increasing pulse duration, the 
maximum charge state as well as the overall signal in- 
tensity grows towards a maximum for an optimal pulse 
width of 800 fs, where atomic ions up to g = 28 can be 
identified. When further increasing the pulse duration, 
both the maximum charge state as well as the overall 
signal degrade. 

The efficient charging for a certain pulse duration was 
in most cases attributed to resonant heating. Another 
mechanism, i.e., enhanced ionization (ENIO), has been 
proposed by (Sicdschlag and Rost, 2003). It relics on the 
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However, because of the large outer Coulomb barriers in 
highly charged clusters, ENIO is considered to be relevant 
primarily for small compounds (A^ ~ 10). Enhanced ion- 
ization due to resonant heating (plasmon-cnhanced ion- 
ization), on the other han d, applies to clusters of any size 
and even to nm-particles ([Doppner et al. ^ 2005 ;_ Reinhard 
and Suraucj, 2001; Saalmann, |200(;; Saalmann and Rost, 



2003; Suraud and Reinhard 



2000) 



conrtept of charge-resonance 
known from diatomic molecules, see Secll 



The emission spectra discussed so far correspond to 
experiments with one single pulse of variable duration. 
Results from a dual-pulse experiment on Ag clusters are 
given in Fig. ^ showing the yield of Ag^'^+ and the maxi- 
mum energy of the emitted electrons as function of pulse 
separation. The ion signal shows a pronounced maximum 
for delays of about 5 ps with an enhancement of more 
than one order of magnitude. This indicates that cluster 
activation and enhanced ionization can be dise ntangled, 
as is also supported by numerical si mulations ( Bornath 



enhanced loni^ion jCni^ij ^/J , ^007^,|b|; |Doppncr et al 



which was 



also considered for multiple ionization of clusters in (Last 
and Jortnci , 1998). Within ENIO, the increased ioniza- 
tion probability occurs for an optimal interatomic dis- 
tance, where the tunnelling barrier between neighboring 
ions and the outer cluster Coulomb barrier of the sys- 
tem are reduced at the same time. The optimal pulse 
duration is thus related to the instant at which the ex- 
panding cluster reaches the optimal interatomic distance. 



2006; Martchcnko et al 



2005 ; Sicdschlag and Rost , 2005). Clear indications that 
a sequence of two pulses may even represent the optimal 
pulse profile for high-q ion production have been found 
by (IZamith et al\ , |00|) for Xbn and by (^ uong et al\. 
200E) for AgAT. ^In both studies genetic feedback algo- 



rithms have been used to optimize the temporal pulse 
structure in order to maximize ion charge states and 
converged towards a pulse profile containing two sub- 
pulses. Dual-pulse excitation further offers a route for 
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targeted control of the cluster dynamics. It has been 
demonstrated on small silver clusters in helium droplets 
that the optimal delay can be controlled by the inten- 
sity of the leading pulse. As was corroborated by semi- 
classical Vlasov calculations, a higher intensity of the 
leading pulse enhances the cluster expansion speed due 
to higher heating and ionization and thus reduces the 
time for which resonant coupling conditions arc estab- 



lished (Doppncr et at, 200E ) 



Coming back to Fig.|36|, also the maximum electron 
energy is analyzed as function of pulse separation. The 
coincidence of high ionization yield and maximal elec- 
tron energy underlines the leading role of collective ex- 
citations in both decay channels. A similar correlation 
bet ween fast electrons an d VUV-radiation was reported 
in (springate et ai, 2003). A common feature is the oc- 
currence of high electron energies. A maximum value of 
375 eV 60 Up was observed with AgAr at moderate in- 
tensity, see Fig. ^ For Xe^y energies in the keV range 
have been reported (|Kumarappan et aL| , 2002; 3hao et al. 



1996; Springate et at, 2003). Further details on the elec- 
tron emission, i.e., angular- and time-resolved signatures 
and underlying acceleration mechanisms arc considered 
in Scc. |VI.B.3| . 



2. Difficulties of explaining high charge states 

Although most of the above trends like higher ioniza- 
tion and energetic particle emission for resonant clus- 
ter excitation can qualitatively be explained, the quan- 



tita1 ive understanding of the emission spectra remains 
a challenge. A still largely debated topic is the ori- 
gin o f the very high atom i c ionization stages from clus- 
ters ( [Fennel et~al\ , |2007b| ; [Heidcnrcich et al\ , ^007| ). In 
order to calculate realistic ion spectra, inner ionization 
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FIG. 36 Comparison of the Ag^°^ yield (diamonds, left 
axis) with the maximum kinetic energy of the emitted elec- 
trons (dots, right axis) following laser excitation of Ag clus- 
ters ((A'^) ^ 2 X 10'') with dual 100 fs laser pulses at intensity 
8 X 10^^ W/cm'^ and 800 nm wavelength. The curv es are guides 



to the eye. Adapted from (Doppner et ai, 2006) 



in the presence of local fields, outer ionization dynam- 
ics, as well as recombination effects have to be taken 
into account consistently. Difficulties arise from at least 
two facts. First, since inner ionization cannot be treated 
fully quantum mechanically for practical reasons, sim- 
pler approximations like ADK-ratcs and atomic impact 
ionization cross-sections have to be used and must be 
corrected correspondingly. Second, recombination pro- 
cesses, even if treated only classically or with effective 
rates, proceed at much longer time scales than the inter- 
action with the pulse and are thus numerically extremely 
time-consuming. However, a few routes towards a more 
realistic description of high charge states by incorporat- 
ing these effects have already been explored. 

To cope with the first problem, inner ionization has 
to be corrected for medium contributions like screening 
or polarization effects in the cluster. This is more or 
less straightforward for tunnel ionization, as the effective 
local field resulting from the momentary distribution of 
charges and the laser is accessible numerically, e.g., from 
MD simulations. Applying an appropriate temporal or 
spatial filtering, the effective field can be used for the 
determination of tunnelling ionization probabilities from 
the ADK-rates (Ammosov et ai, 1986). More involved is 



the treatment of electron-impact ionization which is of- 
ten d escribed by the empirical cross-sections from ( Lotz , 
1967). That requires atomic ionization thresholds which 



are modified by many-particle effects in the cluster, such 
as screening and fields from neighboring ions. One way 
to the determine these shift is the use of statistical ap- 



proaches like Debye screening or ion sphere models (Bor- 



nath et al, 2007a; Gets and Krainov, 2006). These, how- 



ever, assume local thermal equilibrium and neglect the 
details of ionic correlation. A more direct approach re- 
lies on the evaluation of the local field and the resu lting 



shifts dir e ctly fro m a particle based simulation, see (Fcn- 
2007b| ). Irrespective of the particular method. 



nel et al. 



threshold lowering induces substantial enhancement of 
impact ionization when compared to the bare atomic 
cross-sections. An example will be discussed below. 

The second problem concerns the handling of electron- 
ion recombination. Usually it is assumed that only con- 
tinuum electrons produced during the laser pulse con- 
tribute to the final ion spectra and cluster-bound elec- 
trons (quasi-free after the laser pulse) fully recombinc. 
With this assumption, however, the high experimental 
charge states at moderate laser intensities cannot be ex- 
plained. Under experimental conditions this full recom- 
bination of quasi-free electrons is questionable, as, in par- 
ticular, weakly bound electrons may not relax to lower 
ionic levels but can be re-ionized by space-charge field in 
the interaction zone or by ion extraction fields required 



2007b) 



for the time-of- flight analysis (Fennel et al. 

In the latter work it was found that the combined ac- 
tion of both, enhancement of electron impact ionization 
trough threshold-lowering and background-field induced 
frustrated recombination, enhances the maximum ion 
charge states by up to a factor of two, see Fig.p^. While 
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FIG. 37 Calculated maximum charge state of atomic ions 
from Xejv (A*' = 147 — 5083, as indicated) exposed to 250 fs 
laser pulses with peak intensity 4 x 10^''W/cm^. The results 
correspond to different treatments of electron-impact ioniza- 
tion (EII) and electron-ion recombination. Conventional and 
enhanced EII correspond to atomic and local-field corrected 
ionization thresholds, respectively. Total recombination as- 
sumes that cluster-bound electrons recombine with the clos- 
est ion after the laser pulse, while the long-term dynamics 
of quasifree electrons in the presence of a weak static back- 
ground field of 3kV/m is take n into accouiit for f rustrated 
recombination. Adapted from (Fennel el al, 2007b). 



enhanced charging of small clusters is dominated by 
threshold lowering effects, the consideration of the recom- 
bination dynamics becomes increasingly important with 
large clusters. Further contributions such as excitation- 
autoionization or ionization via intermediate states, the 
importance of which is known f o r atomic ele c tron- impact 
ionization ( [Griffin et alj |1984 |Loch et al.[ ^008| ), have 
not been studied in detail yet. 



for A r jv ( Hirokanc et al. , ^004 ; [Kumarappan et al. , 
|200l| ). Fig. |38| displays an example for the directional 
asymmetry for Arjv {N ~ 40000) after excitation with 
(8 X 10^^ W/cm^ ; 100 ps ; 806 nm) pulses, where a polar- 
ization induced shift of the "knee" in the energy dis- 
tribution of about 20% is observed. Similar shifts be- 
tween 15 % and 4 %, dep e nding on pulse duration, were 



found with Xcat ( Li et al. , 2005 ; Springatc et al. , 20001 ) 
and with molecular (N2)jv clusters 



j ( [Krislm amurthy et al 
y|, P006|y 



2004; Mathur and Krislmamurthy 

At least three fundamentally different co ntributions to 
this asymmetry have been described. In (Ishikawa and 



Blenski, 2000), a mechanism is proposed where the ad- 
ditional acceleration is a direct result of the laser field. 
Since the net effect of the laser averages out for ions 
with constant q, rapid charge state oscillations of surface 
atoms were proposed, such that higher effective charge 
states appear during laser half-cycles with outward elec- 
tric field component. That accumulates maximum re- 
pulsion along the laser polarization axis. However, this 
mechanism is unlikely to fully explain the experimental 
asymmetry since the rates for electron-ion recombina- 
tion are very low at the typically high electron tempera- 



tures (Bethe and Salpetcr. 1977) 



The second mechanism is an asymmetric Coulomb ex- 
plosion due to angular depende nt charging of io ns and 
was originally discussed for Cgo (Kou et al, 2000). Near 
the cluster poles, i.e., the regions with surface normal 
parallel to the polarization axis, higher peak electric 
fields from the laser and the cluster field (polarization 
and/or space charge) enhances inner ionization. Thus, 
ions located in this region experience stronger Coulomb 
repulsion. This view of enhanced ion acceleration along 
the polarization axis is supported by numerical simula- 
tions (Fennel et al., 2004; Jungrcuthmayer et al, 2004) 



3. Asymmetric ion and electron emission 

An interesting direction for possible applications of 
clusters is the pulsed generation of energetic ions and 
electrons. The quest for a detailed understanding of the 
acceleration mechanisms is therefore not only driven by 
fundamental interests. The presence of asymmetries in 
angular resolved ion spectra reveals, that the cluster dis- 
integration notably deviates from an isotropic explosion 
process. Further, for excitation with appropriate pulses, 
the electron spectra show strong signatures from field- 
driven acceleration with high directionality. Correspond- 
ing signatures from experimental and theoretical studies 
as well as the main concepts for their explanation are 
reviewed below. 
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FIG. 38 Angular dependence of the ion energy recoil spec- 



a. Angular resolved ion emission: Ion energy spec 
tra exhibit a clear asymmetry, where higher kinetic en- 
ergies appear for the emission along the laser polar- 
ization axis. This w as ffist reported for X cjv (jKuj from ([Kumarappan et al.[ [200 1[) 



marippan et at, 2002; Springatc et al., 2000a) and later 



tra of Ar clusters (A'^=4 x 10 ) exposed to pulses with 
(8x 10^^W/cm2;100ps;806nm). At 0°, the polarization 
of the laser is parallel to the time-of-flight axis. Adapted 
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and the observati on of asymmetric ion charging (Hi- 



rokg jnc et al{ |200J) 



et M ^005| ; [Fennel et d. 



ion acceleration ( 


Breizman 




2004; 


Kumarappan et al, 



200j). In terms of a simple rigid sphere model, clus- 
ter ions and electrons can be described by two homoge- 
nously charged spheres of opposite charge density and 
equal radius. The laser-driven oscillation of the electron 
cloud results in a nonvanishing asymmetric contribution 
to the radial component of the electric field at the cluster 
surface, whereby enhanced repulsion follows for surface 
ions near the cluster poles (Breizman et al, 2005). Even 



for an isotropic ion charge state distribution, this mech- 
anism supports enhanced repulsion for surface ions near 
the cluster poles. This repulsion is particularly strong 
for large-amplitude oscillat ions of the electron cloud at 
resonance (see Sec. VI.B.l[ ). Thus, this model can also 
ex plain the pulse-length-depe ndent asymmetry observed 
in (Kumarappan et al, 2G02| ). 
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FIG. 39 Photoemission spectra from silver clusters 
(A'^ ~ 10'^) exposed to 100 fs laser pulses with peak intensity 
8 X lO^^W/cm^ at 800 nm wavelength: (a) Energy-resolved 
emission parallel (Ej|) and perpendicular (Ex) to the laser 
polarization axis for excitation with a single pulse and dual- 
pulses with optimal temporal delay of At = 1.5 ps. (b) In- 
tegrated signals for three electron energy intervals (as indi- 
cated) and normalized to the maximum obtained for (Ex) as 
function of pulse delay. 



b. Angular resolved electron emission: Compared 
with the ions, the degree of asymmetry is much more pro- 
nounced with electrons. The emission is al i gned t o the 



2002 



Shao 



laser p olarization axis (Kumarappan et al. 
et q|I], 1996| ; springate et al, 2003). This preferential 
ejection is a direct marker for laser-assisted and non- 
thermal emission and turns out to be strongly dependent 
also on the pulse duration. On Xcjv (Kumarappan et al, 
2003a ) find a yield ratio yj|/lj_ ~ 3 for optimal pulse du- 
rations, while almost isotropic emission and less energetic 
electrons are observed for the shortest and most intense 
pulse. The authors related this effect to resonant col- 
lective enhancement of the polarization field. Enhanced 
asymmetry for optimal pulse conditions i s also supported 
by simulations ( Martchenko et al. , 2005 ). 

A pronounced resonance effect has furt her been ob- 
served in a dual-pulse experiment on Agjv ([Fennel et al 



2007a). Two pulses with optimal separation yield simul- 
taneously higher electron energies and stronger asymme- 
try when compared to single-pulse excitation, see Fig[3S|a. 
Comparison of parallel and perpendicular electron yields 
for different energy windows as a function of pulse de- 
lay (Fig.[39|b) shows that the asymmetry increases with 
electron energy. The strongest anisotropy of about 6.5 is 
found for the most energetic electrons, uppermost panel 
in Fig.^9|b. For all chosen energy windows a maximum 
yield is observed for similar delays, supporting the pres- 
ence of plasmon-enhanced electron emission. VUU cal- 
culations on a Nai47, presented in Fig. |4^, show the same 
qualitative behavior. Off-resonance excitation induces 
low-energy electron emission and only a small asymmetry 
(Fig. ^0^ and c) , while resonant dual-pulse excitation re- 
sults in energetic electrons and stronger preference along 
the polarization axis (Fig. ^0[b). Note that both experi- 
ment and calculation show electron energies beyond 60 Up 
along the polarization axis for resonant excitation. 



By trajectory analysis it can be shown that rescat- 
tering of electrons by the cluster potential is crucial for 
the high-energy part of the spectrum. More specifi- 
cally, corresponding electrons gain the major energy frac- 
tion within their final transit through the cluster. The 
importance of rescattering is well-known from atomic 
strong-field ionization. A maximum electron energy of 
10 Up results from backscattering of tunnel-ionized elec- 
trons upon re-encount er with the mother i on at optimal 
laser phase, see, e.g., ( Walker et al , [1996| ). In contrast 



to that, a quasi- linear transit along the laser polariza- 
tion axis with weak deflection turns out to be optimal 
in clusters. The energies can by far exceed the 10 Up 
cut-off from atomic backscattering. Two major effects 
contribute to the energy capture in clusters: (i) accel- 



eration by polarization- fields ( Fennel et al , 2007a), and 
(ii) l aser-field-driven acceleration ( [Saalmann and Ro"st , 
2008). Within process (i), transit electrons travel in 
phase with the dynamic cluster polarization field pro- 
duced from plasmon oscillations. A continuous increase 
of single-particle energy can be accomplished for fully 
matched trajectories. This process of surfacc-plasmon 
assisted re-scattering in clusters (SPARC) supports pref- 
erential ejection of fast electrons along the laser polariza- 
tion axis. It further provides an explanation for strong 
acceleration at the instant of resonant plasmon driv- 
ing due to redistribution of collectively absorbed energy 
to SPARC-electrons. In the simulation run of Fig.|§3, 
the peak amplitude of the polarization field gradient ap- 
proaches 35 GeV/m. This value corresponds to an ef- 
fective intensity 25 times higher than that of the laser 
pulse. Process (ii) results from the laser-driven accelera- 
tion of electrons in a static cluster potential. For a pas- 
sage through a deep global cluster potential, electrons 
acquire high transit velocities. If the velocity and the 
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FIG. 40 Angular resolved electron emission spectra from 
Nai47 exposed to 25 fs laser pulses (800 nm) with peak inten- 
sity 8 X 10^'^W/cm'^, as calculated from semiclassical VUU- 
MD simulations. The data correspond to single-pulse (a) and 
dual-pulse excitations at optimal (b) and a longer nonreso- 
nant delay (c). The emission angle 9 is given with respect 
to the laser polarization a xis. The plot is based on the data 
from (Fennel et al, 2007a) but shown with a more convenient 
intensity scaling. 



polarization axis are parallel and the transit occurs dur- 
ing a beneficial laser half-cycle, electrons can be strongly 
accelerated by the laser field. Also such type of energy 
capture from rescattering, which is most effective with 
a deep cluster potential, produces an alignment of fast 
electrons. Assuming the formation of a particularly deep 
cluster space-charge potential for resonant collective elec- 
tron excitation, this process can result in a plasmon en- 
hancement of the electron kinetic energies as well. A 
detailed analysis and a corresponding sca ling-law for the 
attainable electron energy are given in (3aalmann and 
Ros^, 200S). Besides possible contributions from addi- 
tional many-body effects, the dynamics will contain a 
mixture of the processes (i) and (ii). Nevertheless, mech- 
anism (i) dominates for strong collective motion, e.g., in 
metallic systems at moderate intensity, while (ii) prevails 
with very deep cluster potentials and high laser intensity. 



VII. PERSPECTIVES OF LASER-CLUSTER RESEARCH 

The previous sections have shown that the field of 
laser-irradiated clusters is in an actively developing state. 
Rather than concluding we prefer to discuss a few promis- 



ing future directions. Among those are prospects of 
laser pulse shaping or forthcoming novel light sources. 
Furthermore, complex environments and heterogeneous 
atomic compositions as well as the use of clusters for 
relativistic particle acceleration may open new routes for 
technical applications. Finally we address some prospects 
and challenges of future theory developments. 



A. Laser pulse shaping and control 

One of the intriguing perspectives of light-matter 
coupling pertains to its active manipulation by shap- 
ing th e pulse in amplitude a nd ph ase 



et al 



2001a; Brixner and Gerber, 2003). With molecules 



this approach follows the suggestion of (Judson and 
1992| ), in which a computer-controlled pulse 



shap er is used in combination with a lear n ing alg orithm, 
see (Baunicrt et at, 1997; Brixner et al, 2001b| ), in or- 



der to achieve, e.g., a selective molecular reaction. The 
quantum-mechanical processes can be controlled with the 
direct feedback from the experiment in an automated 
fashion, without requiring any model for the system re- 
sponse. This electron wave-packet engineering has be- 
come a power f ul too l to realize the concept of femtochem- 
istry (Zcwail, 1980 ). A more recent technological devel- 
opment further increases the possibilities and prospects 
of quantum control. With the technique of femtosec ond 
polarization pulse shaping ( Brixner et al , 2002 , 2004), it 
is now possible to vary intensity, instantencous frequency, 
and light polarization (i.e., the degree of cUipticity as well 
as the orientation of the principal axes) as functions of 
time within a single femtosecond laser pulse. Thus, a full 
temporal and spatial control is at reach. 

For intense laser-cluster interactions, shaping the pulse 
in amplitude and phase can be a fascinating tool to selec- 
tively steer the dynamics of charging, particle-, or pho- 
ton emission. Basic findings along this line are the con- 
trol of the Coulomb explosion by varying the laser pulse 
length as well as the time delay in the dual-pulse exper- 
iments as outlined in Sec. VI.B.l. For example. Fig. ^ 
has shown the dramatic effect of the adjustment of two 
femtosecond laser pulses on the charging efficiency and 
the energy of emitted electrons. Adaptive femtosecond 
contro l was demonst r ated on the Coulomb explosion of 
Xcjv ( Zamith et al , 2004). Here the signal of highly 
charged Xe''+ could be optimized with the help of a 
simple genetic algorithm, applied to an initially Fourier 
transform limited pulse with 100 fs duration and 230 /iJ 
energy. The procedure resulted into a sequence of two 
120 fs pulses with similar amplitude and separated in 
time by abo ut SOOfs, like in the optimized dual-pulse 
experiments (Doppner et al, 2005). It is very interest- 
ing to note that this two-pulse optimum has been worked 
out by the algorithm starting from a 80 parameter un- 
biased configuration. Corresponding simulations within 
a semiclassical molecular dynamics approach predicted 
that for selected combinations of cluster size, laser inten- 
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sity, and wavelength, th e ionization may b e optim ized by 



a three-pulse sequence ( [Vlartchcnko et al. , 2005 ). In an- 



other closed-loop optimal control experiment on rare gas 
clusters, pulse-shaping has s hown a significant p otential 
for x-ray yield enhancement (Moore et al, 20G5| ). 

Whereas the optimal-control studies on clusters were 
limited to an optimization of the pulse amplitude so 
far, the simultaneous variation of the pulse phase is still 
an exciting challenge. First results of such a fully un- 
biased adaptive fs experiment have demonstrated the 
controlled adjustment of charge state distributions from 
the Coulomb explosion of Ag^r embedded in Helium 

[200S| ). 



droplets (Truong et al. 



In this study the opti- 
mization of the Ag'^+ charge spectrum converged to a 
pulse structure with a weaker pre-pulse and a stronger, 
negatively chirped main pulse. However, we are far off 
a fu ll theoretical understanding of the complex dynam- 



ics ( [IrivcQ by pulses shaped in amplitude and phase. In 
the future, if sufficient mass-selected cluster intensity can 
be prepared, single ionization states and narrow-banded 
high-energy radiation might be realized. 



B. Towards VUV-, XUV, and soft x-ray pulses 

The nature of the laser cluster coupling fundamentally 
changes when going from the IR regime towards exci- 
tation with VUV-, XUV-, ore even x-ray pulses. This 
concerns ionization processes as well as the mechanisms 
of energy absorption. For excitation with IR pulses, field- 
driven ionization plays a crucial role for the nanoplasma 
generation, e.g., in rare gas systems. The subsequent 
energy capture, which eventually removes electrons from 
the cluster, is of plasma nature and can be strongly en- 
hanced through resonant collective excitations. Because 
of extensive plasma heating and resulting further ioniza- 
tion high charge states can arise with IR pulses. 

When going below about 100 nm wavelength, a value 
which was used in the first VUV experiments on rare gas 
clusters, photoionization becomes the dominant charg- 
ing mechanism for inner ionization. Concerning the en- 
ergy absorption, collective effects can be disregarded as 
the required critical density cannot be reached and pure 
IBS heating prevails. In fact, the observation of surpris- 
ingly high energy ca pture in the first VUV-experiments 
on clusters, sec Sec. VLA.2 , have sparked a remarkable 
progress in the understa nding of heating- and ionization 



effects in dense targets (Bauer 



2007b|; Jungrcuthmayer et al. 



2005 



2004b; peorgcscu et al. 



et al. 



20061: ISaalmann et alj |2006|; [Santra and Greene| , |2003 



Siedschlag and Rost|, |2004|; ^iaja et al\, ^007\) 



\ Micn further increasing the laser frequency, IBS heat- 
ing becomes more and more suppressed, cf. Eq. (|^), 
so that photoexcitation of tighly bound electrons be- 
gins to become even the leading energy capture pro- 
cess. Signatures of this transition have recently been ob- 
serv ed on Arjy in intense femtosecond XUV FEL pulses 
at A | = 32 nm {hcji^s 38 eV) ([Bostcdt et al[ |200S|). By 



comparing the experimental photoelectron spectra with 
complementary Monte-Carlo simulations, the following 
behavior was found. The cluster ionization first proceeds 
as a multistep process of direct single-photon absorption 
events. Electrons are released from the cluster directly 
without prior inner ionization and the space charge build- 
up results in an energy downshift for subsequent ion- 
ization steps. This shift leads to a highly nonthermal 
electron energy distribution. At a certain degree of ion- 
ization, the cluster potential frustrates further electron 
release, leading to the formation of a nanoplasma only 
beyond a certain threshold intensity. Even at higher in- 
tensity no strong impact of IBS heating was found. These 
findin gs are in agreement with corresponding MD re- 
sults ( Arbeiter and Fennc] , 2009) and calcu l ation s based 

on kinetic transport equations (Ziaja et al., 2009|}^ 

Using int ense soft x-ray pulses at A = 13 nm (prioener 



et al. , 2008 ) found highly efficient charging of Xejv with 



ions uplC?+ , which can be ascribed to the large absorp- 
tion cross section of the giant atomic Xe A-d resonance. 
By surrounding Xe^r with an additional argon layer it 
was further shown, that charge recombination dynamics 
can be studied in the well controllable core-shell system. 

Another interesting issue concerns the time-resolved 
monitoring of the cluster excitation and the subse- 
quent Coulomb explosion by combining different types 
of pulses. For example, the ionization of rare gas clus- 
ters imy_be_dTiven_by VUV radiation, as in the case 
of (Wabnitz et al.., 2002), whereas a subsequent IR pulse 
probes the collective electron response of the priorily met- 
allized system (Siedschlag and Rost, 2005). A combi- 
nation of VUV and XUV pulses was proposed to mon- 
itor the time-depend e nt ion ization stages in small clus- 
ters ( Georgescu et al. , 2007a ) . Another scheme uses x-ray 



radiation for Thomson scattering on exploding clusters or 
droplet s, which have be en initially excited by strong IR 
pulses ( Holl et al , 2007). By this, a fundamental under- 
standing can be gained on highly nonstationary strongly 
coupled plasmas and their transition from degenerate to 
classical systems. The advent of the X-FEL will open di- 
rect access to the temporal development of such complex 
systems. 



C. Clusters In an environment 

Embedding clusters into an environment or deposit- 
ing them at surfaces modifies their optical responses, 
sec ( pCreibig and VoUmci , 1995| ). A major branch of 
present-days cluster research comprises system s in con- 



tact with so l id su r faces, for a review see, e.g . ,{ Meiwes- 



Broer, 2000, 2006; Meiwes-Broer and Bcrndt, 2007). An 



extremely rich scenery unfolds when considering the spe- 
cific effects emerging from the interaction of a cluster 
with an environment. One finds, e.g., only small shifts 
for the Mie plasmon res onances of meta l clust e rs em- 
bedded in in ert matter ( Diederich et al. , ^002 ; Fehrei 



et al, 2007b) and larger ones for contact with conduct- 
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FIG. 41 Time evolution of the radius mean square 
(rms) of free Nag (bottom), Nag embedded in Ar434 



(middle), and 
(2.4 lO^nV/c 



pure 
33 fi 



Ar447 (top), 
OOOmn). 



arter n'radiation with 



/ cm , ib, u jumn j . Calc ulaLiuns liave been 
performed using TDDFT for Na valence electrons an d MD 
for Na+ ions and Ar atoms. From (Fehrer et al, 2007a), with 



kind permission of The European Physical Journal (EPJ) 



ing material ( Pinchuk et al , ^004 ). Details of the exci- 
tation spectrum, however, are rather sensitive to the in- 
terface. This affects, e.g., the spectral fragmentation and 
the p lasmon damping , for e xperimental a s sessm ent see, 
e.g., (Hendrich et al, 2003; Zicglcr et al, 2004). Large 
effects from the environment are to be expected in the 
reaction dynamics at high excitations. 

A theoretical example is shown in Fig. It compares 
three test cases, Nas as a small metal cluster, the Nas 
embedded in Ar434 (a large rare gas cluster as model for 



a m atrix), and pure Ar44;, all three cxpoocd to the oame 



lase r puloo. — The laoor puloo loado to a charge otato 3 I 
of Nag. In the free case (bottom panel), this induces 
a Coulomb explosion. The situation is quite different 
for Nag in Ar434 matrix. The metal cluster is again 
highly excited and starts to explode. But the explosion is 
stopped by the Ar atoms which efficiently absorb the ex- 
citation energy of the system (lower middle panel). The 
Ar matrix is perturbed and exhibits monopole oscilla- 
tions, but of much smaller amplitude than Nag (upper 
middle panel). The upper panel of Fig. ^ finally shows 
the case of a pure Ar447. Under the same laser conditions, 
one can see that the Ar447 remains essentially unper- 
turbed, showing no electron emission and only extremely 



weak breathing oscillations. Obviously, the Nag acts here 
as a chromophore, absorbing energy from the laser pulse 
and transferring it to the environment. The example 
shows that the combination of two materials changes the 
reaction dynamics of either system dramatically. One can 
easily imagine that putting clusters in contact with var- 
ious substrates produces a world of interesting scientific 
questions and offers technical applications in the field of 
nanotechnology. We mention in the following briefly a 
few aspects to give an idea of the enormous possibilities, 
concentrating on optical properties. 

When depositing Auat on a semiconductor surface the 
change of optical cluster propertie s can be explo i ted t o 
producing enhanced photocurrent (Schaadt et al, 2005). 
There are promising applications, e.g., in medicine where 
the frequency selective optical coupling of organically 
coated metal clusters attached to biological t issue may 
by us e d for d iagnosis ( Bruchcz Jr et al ^ 1998|; [Dubertret 



et al\ , |2002| ; [Mayer et al\ , |2001| ; pimberg et al[ |2007l ) 



or, in the case o f stronger laser fi e lds, fo r localized heat 
ing in therapy (Khlebtsov et al, 200(;). The field am- 
plification effect is of interest in many other materials 
and applications, see, e.g., the study of localized melt- 
ing for the generic combination of Au clusters embedded 
in ice (Richardson et al, 2006). The strong coupling to 
light may be used for more than just heating. Ensembles 
of size and shape selected c lusters on a surface are pro- 
duced by "s hape burning" ( Ouacha et al , 2005| ; Wenzel 
et al, 1999| ). A dedicated modification o f the shape for 



embedded Ag clusters is demonstrated in (Dahmen et al 



2006; Perner et al, 2000). Time scales and mechanisms 
of energy transport are thus an issue of high interest. 
Theoretical analysis has yet to cope with the great vari- 
ety of material combinations. For an example using the 
generic t est system of metal cluste r in a rare gas ma- 
trix, see ( [Fehrer et al] , [2007a[ , ^008[ ). A thorough study 
of surface-deposited cluster subject to strong laser pulses 
still is a matter of future studies. 



D. Relativistic particle acceleration with clusters 

Strong laser fields impinging on clusters can drive in- 
teresting electron dynamics. For an example from the 
mode r ate in tensity domain (^ lO^^W/cm^), (Fennel 



et al 



^007a ) describe a cascade- like acceleration mecha- 
nism based on resonant field amplification in individual 
clusters, see Sec JVI.B.3| . In the regime of 10^^ W/cm^ 



electron energi es from keV u p to some hundreds o f kcV 



al], ]2002a[ ; ^hao et alj [1996 



are reported (Chen et 

Bpringate et al , 2003[ ), emitted in transverse direction to 
the laser propagation axis. Beyond a few tens of keV the 
emission is most likely due to macroscopic plasma wave- 
breaking effects in a very dense cluster beam, as is further 
supported by a pronounced forward peak in the emis- 
sion (Chen et al, 2002b). Moreover, there are few exam- 
ples close to or in the relativistic regime (lO^^W/cm^). 
From studies on bulk and dense atomic gases it is known 
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that charged particles can be accelerated by the plasma 
Wakefield to large k metic energies, for a detailed t h eoret 
ical discussion see (Pukhov and Meyer-ter-Vchn 



200 



and for a recent experimental example ( Karsch et al. 
2007). There exist realistic plans to employ the effect 
to buil d fairly inexpensive laser-driven table-top acceler- 
ators ( priincr et al. , 2007). Indications for special rcl- 
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FIG^iS — ElycLiuu kliiyLlc eiiyigy dlbLiibuLluii lysulLliig fiuin 
irradiation of large Ar particles (micron sized diameter) in Ar 
gas environment by an intense laser (3.5x lO'^^W/cm^). The 
solid, nir-ifp shows the detprtinn threshnld Strflight lines in- 



and dotted for T = 2.8 MeV. Adapted from (Fukuda et al. 
2007), with permission from Elsevier. 



ativistic electron acceleration mechanisms with clusters 
have been reported in a stud y on large Ar p articl es in a 



20071) . The 



low-density background gas ( Fukuda et al. 
example in Fig 



2^ displays the achieved electron kinetic 
energies. On the basis of simulation results, the two tem- 
peratures have been associated with two different gener- 
ating mechanisms. The lower-energy electrons stem from 
acceleration in a distorted wakefield. In contrast to that, 
the high-energy electrons are removed from the clusters 
with already relativistic energies and then further accel- 
erated by the laser pulse directly. The kinetic energies 
observed here are still far below what can emerge from 
bulk plasma. However, whereas optimum conditions, ad- 
vantages, and disadvantages have yet to be worked out, 
the example proves the feasibility of relativistic particle 
acceleration with clusters. Moreover, the use of clusters 
as dense electron containers for free-space electron accel- 
eration, e.g., with radially polarized laser beams, might 
be promising for generating ultrashort electron bunches 
with durations down to the attosccond domain at up to 
GeV e nergy (Karmakar and Pukhov, 2007; Varin and 
Pich4 pOM). 



E. Challenges for theory 

The theoretical description of laser-cluster dynamics 
requires to cover very different scales of length, time, and 
energy. This difficulty usually hampers a fully micro- 
scopic treatment of all degrees of freedom. Fortunately, 



the nature of the response is to great extend determined 
by the type of excitation or the size and structure of 
the tar get, at least at the different scales. We have seen 
in Sec. Ill that there exists a bunch of theoretical ap- 
proaches, ranging from fully microscopic ones to macro- 
scopic ones, which are applicable within certain windows 
of size and energy. Their limitations result from both, 
formal constraints, e.g., due to the level on which corre- 
lations and quantum effects are resolved, as well as from 
practical ones like the numerical expense. To promote 
the development of more elaborate methods and schemes 
with wider ranges of applicability we see at least two 
promising directions. 

The first and most straightforward path is the formal 
improvement of particular methods along with the rapid 
development of high-performance computers. For in- 
stance, the impressive growth of numerical power allows 
the application of fully correlated quantum approaches 
to systems with several electrons, e.g., with the efficient 
handl i ng of few-body wavefunction s by MCTDH 



et al\ |00g) or MCTDHF methods ( |Caillat et al. 




This opens a route to explore truly correlated electron 
dynamics including continuum and intermediate excited 



dica le fits to thermal distributions, dashed for T = 18.8 MeV asyan et 



states, for an exam ple on a molecular system see ( 3uki- 
l2009| ). A mental challenge con- 



cerns the inclusion of dynamical correlations in mean- 
field quantum theories like TDLDA in the sense of a 
quantum counterpart to the semiclassical description 
within VUU. On the classical level, efficient numerical 
schemes and large-scale parallelization promises the fea- 
sibility of up to giga-particle simulations. 

A second frontier concerns the connection of different 
treatments in terms of multi-level or multi-scale methods. 
A well-known example for biological and chemical appli- 
cations are mixed quant um mechanics / molecular d ynam- 
ics approaches, see e.g. (Bakowies and Thiol, 1996| ). Also 
for clusters the combination of different levels has turned 
out to be very successful, e.g., within MD and hydro- 
dynamic schemes for strong-field excitations, where the 
quantum nature of inner ionization is taken into account 
via effective rates and cross-sections. The connection of 
different treatments, however, requires interfaces, the val- 
idation of which is a great challenge. Firm links between 
the approaches and reliable interfaces, e.g., within over- 
lapping zones like those indicated in Fig. ra, are therefore 
highly desirable and have far-reaching implications. One 
example could be the connection of an explicit atomic- 
scale quantum treatment of inner ionization with a more 
coarse-grained semiclassical or even classical treatment 
of quasifree and continuum electrons. This would be of 
high interest for strong-field laser-cluster interactions in 
a wide range of laser frequencies, i.e., from the IR up to 
the x-ray domain. Another challenging aspect are strong- 
field excitations of larger clusters and particles in the IR 
range, where propagation effects of the light field cannot 
longer be neglected. Here a combination of molecular dy- 
namics techniques for evaluating the short-range part of 
the interactions combined with electromagnetic particle- 
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in-ccU concepts for describing the long-range component 
of the Coulomb- and radiation fields might be promising. 
Last but not least, such neighboring approaches could 
also be combined in a sequential way, e.g., to resolve the 
laser excitation microscopically, whereas the long-term 
behavior is described with a less expensive scheme. 

Along these lines the field of laser-cluster dynamics 
will certainly be inspired by forthcoming developments 
in other branches like atomic, molecular, and plasma 
physics and vice versa. 
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